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I.  STATEMENT  OF  WORK 

Conduct  a  theoretical  research  program  to  develop  quantum  mechanical 
methods  for  studying  nonadiabatic  effects  in  three-dimensional  atom-molecule 
interactive  collisions. 

II.  TECHNICAL  DISCUSSION 

The  traditional  role  of  theory  in  chemistry  has  been  to  develop  models 
that  help  to  correlate  large  amounts  of  experimental  data.  However,  because 
of  the  difficulty  in  obtaining  atate-to-state  data  on  many  systems,  there  is 
growing  interest  in  using  the  computational  methods  of  ab  initio  theoretical 
chemistry  to  predict  kinetic  data.  Due  to  the  increase  in  sophistication  and 
accuracy  of  the  computational  methods  available  for  calculating  potential 
surfaces  and  obtaining  dynamical  information,  theoretical  studies  specifically 
designed  to  predict  kinetic  data  are  appearing.^ 

Many  areas  of  military  technology  are  affected  by  the  new  developments  in 
both  experimental  and  theoretical  kinetics.  The  collisional  mechanisms  of 
importance  in  chemical  lasers  include  the  conversion  of  vibrational  energy 
into  rotational  energy  (V-RT),  conversion  of  rotational  energy  into 
translational  energy  (R-RT),  exchange  of  vibrational  energy  (V-V), 
intramolecular  conversion  of  rotational  energy  into  vibrational  energy  (R-V), 


translational  to  internal  energy  conversion  (T-RV) ,  and  the  conversion  of 
electronic  energy  to  vibrational,  rotational,  and  translational  energy  (E- 
VRT).  Also  of  importance  is  the  transfer  of  electronic  energy  from  one  species 
to  another.  For  example,  the  transfer  of  electronic  energy  from  02(*A)  to 
atomic  and  molecular  iodine  is  important  in  iodine  lasers.  In  interhalogen 
lasers  such  as  IF  the  production  of  electronically  excited  IF  from  F,  I^,  and 
02(*A)  is  poorly  understood. 

In  the  area  of  plume  technology,  an  important  problem  is  the  prediction  of 
the  radiant  intensity  of  rocket  exhaust  plumes  in  the  upper  atmosphere.  These 
plume  signatures  are  driven  by  collisional  excitation  of  exhaust  gas  molecules 
by  ambient  atmospheric  species,  and  standard  models  that  predict  plume 
radiation  require  state-to-state  cross  sections  for  collisional  excitation  as 
input.  The  necessary  cross-sections  potentially  involve  all  of  the  energetic 
processes  which  are  also  of  importance  in  chemical  lasers.  Theoretical 
methods  provide  a  cost  effective  and  timely  means  of  obtaining  data  of  this 
type.  For  systems  which  are  inaccessible  to  current  experimental  techniques, 
theoretical  methods  provide  a  complementary  tool  to  the  experimental 
determination  of  kinetic  data. 

The  dynamical  method  most  widely  used  for  studies  of  complex  systems  is 

3 

the  quasiclassical  trajectory  method.  However,  many  phenomena  cannot  be 
accurately  described  by  a  classical  approach,  and  quantal  methods  are 
required.  Present  quantum  theoretical  methods  for  obtaining  transition 
probabilities  are  difficult  to  apply  to  problems  such  as  the  vibrational 
excitation  of  HF  (or  even  Hj)  because  of  the  large  number  of  internal  states 
that  must  be  included  to  converge  the  calculations.^  There  have  been  several 
exact  close-coupling  calculations  reported  for  inelastic  collisions  involving 
light  three-atom  systems  at  relatively  low  energies,  mostly  involving  rigid- 
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rotor  models,  but  ex.  ct  close-coupling  studies  of  bimolecular  reactions  have 
been  limited  to  H+Hg.  Because  of  the  importance  of  obtaining  quantum 
mechanical  results  for  collision  processes,  approximate  decoupling  methods 
have  been  developed. Besides  greatly  enlarging  the  scope  of  inelastic 
quantum  mechanical  calculations  to  include  vibrational  excitation  processes, 
implementation  of  decoupling  methods  is  allowing  three  dimensional  reactive 
scattering  calculations  for  systems  more  complex  and  of  greater  interest  than 

h+h2.7 

One  area  of  current  research  that  has  received  considerable  attention  is 
the  development  of  methods  for  including  electronic  degrees  of  freedom  in 

O 

scattering  calculations.  The  usual  difficulties  of  the  coupled  channel 
approach  are  magnified  not  only  because  of  the  many  s»re  states  introduced, 
but  due  to  complicated  angular  momentum  coupling  schemes  not  present  in 
studies  on  a  single  adiabatic  electronic  surface.  Thus  there  are  few 
instances  where  actual  applications  to  three  dimensional  atom-molecule  systems 
have  been  made.  Progress  in  this  area  will  involve  improvements  in  techniques 
for  treating  the  large  number  of  equations  involved,  particularly  for  reactive 
systems,  where  the  further  complications  of  difficult  coordinate  systems 
appear . 

This  report  summarises  the  progress  at  Chemical  Dynamics  to  extend  the 
quantum  mechanical  approach  to  molecular  dynamics  to  systems  of  greater 
complexity  than  previously  possible.  The  emphasis  is  on  obtaining  detailed 
infonsation  on  vibrational  and  rotational  nonadiabatic  processes  in  atom- 
molecule  reactions,  and  in  developing  methods  for  including  electronic  degrees 
of  freedom  in  three-dimensional  quantum  scattering  theory. 


III.  RESEARCH  OBJECTIVES 


The  overall  objectives  of  this  research  program  are  as  follows: 


•  Extend  current  technology  to  treat  the  nonreactive,  electronically 
nonadiabatic  collisions  of  an  atom  with  a  diatomic  molecule.  This 
involves  the  development  of  a  code  for  incorporating  the 
appropriate  electronic  basis  sets  and  decoupling  schemes. 


#  Investigate  decoupling  approximations  suitable  for  studying 

vibronic  transitions.  Such  processes  are  very  important  in 

understanding  the  physical  processes  occurring  in  many  areas  of 
technology. 


•  Develop  strategies  for  incorporating  electronic  angular  momentum  in 
molecular  collisions. 


•  Conduct  computational  applications  to  the  systems  O+Hj  and  Na+H, , 
and  other  systems  appropriate  for  the  development  of  tne 
computational  methods.  The  calculations  will  involve  both 
nonreactive  and  reactive  studies. 


•  Continue  the  development  of  quantum  scattering  codes  to  enable  more 
systems  of  chemical  relevance  to  be  amenable  to  study  by  quantum 
mechanical  means.  This  development  is  an  important  component  of 
research  at  Chemical  Dynamics,  where  a  major  goal  is  the 
development  of  state-of-the-art  methods  capable  of  providing  the 
best  theoretical  data  for  problems  of  current  interest  in  chemical 
kinetics. 


In  addition  to  progress  in  computational  techniques,  the  goals  of  this 
program  involve  furthering  our  understanding  of  fundamental  processes  in 
molecular  collisions,  such  as  resonance  phenomena  and  energy  pathways.  Recent 
developments  in  quantum  chemistry  permit  the  calculation  of  accurate  potential 
energy  surfaces  for  both  ground  and  excited  molecular  states,  and  provide  new 
justification  for  the  search  for  more  powerful  methods  for  studying  dynamical 


processes 
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IV.  RESEARCH  ACCOMPLISHMENTS 


During  this  eighteen  month  research  program,  there  has  been  progress  in 
several  areas.  These  include  the  development  of  a  computational  method  for 
analyzing  three-dimensional  reactive  scattering  wavefunctions  as  a  function  of 
reaction  path  and  for  generating  wavefunctions  of  sufficient  accuracy  to  permit 
the  computation  of  angular  distributions  for  direct  comparison  with  experiment. 
Applications  to  the  F+Hj  reaction  are  providing  new  insight  into  energy  disposal 
during  exothermic  reactions  and  into  the  nature  of  resonance  phenomena  (see 
Appendix  A).  The  quantal  angular  distributions  show  features  that  appear  in 
experimental  data,  and  that  are  not  predicted  by  classical  theoretical  methods 
(see  Appendix  B).  Preliminary  close-coupling  probabilities  have  been  obtained 
for  this  reaction  for  J>0.  These  results  are  presently  not  sufficiently 
converged,  even  with  the  180  channels  so  far  included,  but  there  is  encouraging 
agreement  between  these  results  and  the  J^-conserving  probabilities  that  have 
formed  the  basis  for  our  previous  studies  of  this  system.  We  should  point  out 
that  these  results,  obtained  under  AFOSR  sponsorship,  are  by  far  the  most 
sophisticated  quantum  mechanical  studies  of  chemical  reactions  to  date. 

Previous  work  on  this  contract  included  the  development  of  techniques  for 

treating  electronic  transitions  in  molecular  collisions  using  both  adiabatic  and 

diabetic  representations.  These  new  methods  have  been  applied  to  the  problem  of 
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the  quenching  of  the  resonance  state  of  potassium  by  hydrogen  and  muonium.  ’ 
Muonium  (Mu)  is  an  electron-muon  pair  which  behaves  as  a  light  isotope  of 
hydrogen.  Because  of  the  extremely  light  mass  of  muonium  (one-ninth  that  of 
hydrogen)  the  validity  of  the  Born-Oppenheimer  adiabatic  separation  of  nuclear 
and  electronic  motion  has  been  questioned. ^  The  comparison  of  cross  sections 
for  the  quenching  of  K  by  H  and  Mu  has  given  insight  into  the  applicability  of 
the  Born-Oppenheimer  approximation  in  systems  containing  Mu. 


Work  is  currently  in  progress  to  examine  the  rich  resonance  features  in 
this  relatively  simple  atom-atom  collision.  Plots  of  the  cross  sections  versus 
translational  energy  for  the  K  +  Hu  collision  show  sharp  peaks  which  are 
attributed  to  shape  type  resonances  at  low  translational  energies.  This 
characterization  of  the  low  energy  resonances  has  been  confirmed  by  showing  that 
the  energies  at  which  the  resonances  occur  correspond  to  quasibound  states  of 
the  effective  potential  for  the  first  excited  electronic  state  of  KMu.  Similar 
resonance  features  are  expected  for  the  KH  system;  however,  these  resonances  are 
expected  to  be  much  more  narrow  because  of  the  heavier  mass  of  hydrogen. 

Resonances  provide  a  very  detailed  probe  of  the  potential  energy  surface. 

« 

Comparisons  of  calculations  of  these  types  with  spectroscopic  observations  of 
the  resonance  states  of  the  "bound"  molecule  can  provided  information  about  the 
potential  curves  in  regions  previously  undefined  by  tools  such  as  RKR  analysis. 

An  improved  Diatomics- in-Molecules  (DIM)  representation  for  the  lowest 
electronic  surfaces  of  Ha+Hj  has  been  obtained  for  use  in  scattering 
calculations.**’*^  These  surfaces  are  being  used  in  our  current  dynamical 
studies  of  this  system  (see  Appendix  C). 

Significant  improvements  have  been  made  in  the  iterative-variational 
method  for  heavy  particle  scattering  (see  Appendix  D).*^  These  include  the  use 
of  finite-element  techniques  for  greatly  reducing  the  amount  of  I/O  required  by 
the  previous  implementation,  and  the  development  of  a  scheme  for  allowing  for 
chemical  reaction  through  an  efficient  implementation  of  Miller's  exchange- 
kernal  formalism. *  ^  ’*^  To  date,  development  of  an  inelastic  code  has  been 
started,  with  the  implementation  of  the  finite  element  method  complete. 
Implementation  of  an  initial  iterative  scheme  has  been  completed,  a'd  testing  on 
the  problem  of  HF-HF  self-relaxation  begun  but  not  fini~.*».  ,  is  development 
was  suspended  due  to  lack  of  funds.  However,  we  feel  that  this  approach  holds 


promise  for  significantly  expanding  the  number  of  channels  that  can  be  included 
in  quantum  scattering  calculations  and  warrants  further  development.  It  is,  in 
our  opinion,  the  approach  most  likely  to  succeed  in  overcoming  the  limitation  of 
coupled-channel  methods  imposed  by  the  mass-dependence  of  many  molecular 
systems,  particularly  those  with  small  skew-angles  in  standard  reaction 
coordinate  theories.  The  current  reaction  coordinate  methods  are  intractable 
for  molecular  systems  involving  masses  other  than  hydrogen,  except  in  a  few 
cases  where  decoupling  approximations  work  to  an  acceptable  degree.  This 
computational  method,  installed  on  one  of  the  new  generation  of  supercomputers, 
could  open  new  horizons  for  accurate  computation.  With  their  very  large 
memories  and,  more  importantly  for  the  present  application,  very  rapid  and  large 
mass  storage,  these  computers  will  possess  the  capacity  to  apply  such  algorithms 
to  scattering  problems  on  the  order  of  thousands  of  channels.  This  will  enable 
the  accurate  methods  of  quantum  scattering  theory  to  be  applied  to  situations 
that  are  presently  tractable  only  to  classical  methods. 
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V.  RESEARCH  IN  PROGRESS 

During  t^’s  program,  computational  applications  to  several  important 
collision  collision  problems,  including  studies  of  nonadiabatic  processes  in 
Na+Hj  collisions  have  been  initiated.  These  calculations  incorporate  the 
Infinite-Order  Sudden  (IOS)  approximation  and  a  new  approach  to  including 
electronic  coupling.  This  method  allows  the  economical  study  of  electronic-to- 
vibrational  energy  processes  by  eliminating  the  explicit  treatment  of  rotational 
degrees  of  freedom. 

Reactive  calculations  on  0(  P)*!^  are  in  progress,  and,  when  sufficient 
computer  time  is  available,  converged  cross  sections  will  be  obtained.  This 
system  is  of  interest  because  of  the  previous  classical  trajectory  and 
transition  state  theory  treatments,  the  availability  of  a  reasonable  ab-initio 
surface,  and  the  existence  of  experimental  rate  data.  The  present  three- 
dimensional  treatment  is  tractable  because  of  the  relatively  small 
endothermic ity  of  the  system  and  the  presence  of  hydrogen  atoms  in  both  reagent 
and  product  molecules. 

Work  is  also  continuing  on  the  problem  of  obtaining  converged  close¬ 
coupling  probabilities  for  the  F+H^  reaction.  This  is  an  important  problem,  for 
these  probabilities  are  necessary  to  fully  evaluate  the  reliability  of  the 
decoupling  methods  that  are  essential  for  obtaining  complete  cross  sections.  As 
of  the  present  time,  we  have  made  close- coup ling  calculations  for  J  =  2  that 
involve  approximately  180  coupled  channels.  The  basis  of  rotor  functions  has 
been  varied  among  the  lower  vibrational  states  to  study  convergence.  The 
results  compare  reasonably  well  with  the  converged  J  -conserving  calculations  at 
low  energies.  The  latter  include  up  to  140  channels.  The  total  probability  for 
the  process  F+HjCv^O , j*0)  ->  H+FH(v*2,e j)  is  0.39  in  the  .^-conserving 
calculation  compared  with  0.43  in  the  close-coupling  calculation.  The  close- 
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coupling  result  is  still  unconverged.  We  are  presently  in  the  process  of 
rewriting  the  scattering  code,  using  a  new  partitioning  of  the  Hamiltonian,  to 
allow  many  more  channels  to  be  included  so  that  the  close-coupling  calculations 
can  be  converged  at  least  for  small  J. 
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VIII.  INTERACTIONS 
A.  MEETINGS 


The  following  meetings  were  attended,  providing  valuable  interaction  with 
other  scientists  and  DoD  personnel: 

Austin  Conference  on  Theoretical  Chemistry,  March  1981  (M.J.  Redmon  and 
B.C.  Garrett) 


12th  ICPEAC ,  Gatlinburg,  July  1981  (B.C.  Garrett) 


Gordon  Conference  on  Atomic  and  Molecular  Collisions,  Plymouth,  N.H. , 
July  1981  (B.C.  Garrett). 

AFOSR  Dynamics  Contractors  Meeting,  Albuquerque,  October  1981  (M.J. 

Redmon) 

DoD  Plume  Visibility  Workshop,  Huntsville,  November  1981  (M.J.  Redmon) 

AFRPL/AFSD  Plume  Technology  Workshop,  Los  Angeles,  February  1982  (M.J. 
Redmon) 

American  Chemical  Society  Annual  Meeting,  Las  Vegas,  Nevada,  March  1982 
(L.D.  Thomas) 


Symposium  on  Theoretical  Aspects  of  Gas  Phase  Kinetics,  University  of 
Reading,  England,  March  1982  (B.C.  Garrett) 


Faraday  Disussions  on  Van  der  Waals  Molecules,  Oxford,  England,  April 
1982  (B.C.  Garrett) 


Workshop  on  the  fundamentals  of  initiation  of  chemical  decomposition  of 
energetic  materials,  Chester,  Maryland,  May  1982  (B.C.  Garrett) 

American  Chemical  Society  Meeting,  Kansas  City,  Missouri,  September  1982 
(B.C.  Garrett) 
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B.  SUPPORT  OF  OTHER  HOD  PROGRAMS 

The  techniques  investigated  under  the  current  program  have  application  to  many 
areas  of  DoD  interest  involving  gas-phase  chemistry.  He  have  recently 
received  support  form  the  Rocket  Propulsion  Laboratory  at  Edvards  Air  Force 
Base  to  study  collisional  excitation  of  various  high-altitude  exhaust  gas 
species  by  atomic  oxygen.  The  quantum  scattering  calculations  will  be 
performed  with  methods  developed  under  the  current  program.  The  computer 
codes  developed  under  the  current  contract  are  to  our  knowledge  the  only 
existent  codes  capable  of  performing  these  calculations. 

He  have  also  submitted  a  proposal  to  perform  a  theoretical  study  of 
rotational  relaxation  in  HF-HF  collisions,  a  subject  of  enormous  practical 
importance  to  the  Air  Force  HF  laser  program.  This  study  would  involve 
application  of  the  iterative-variational  method  in  a  close-coupling  mode  for 
the  loir’ lying  v-1  states  decaying  to  the  v*>0  manifold,  with  decoupling  applied 
to  higher  states.  There  is  excellent  experimental  data  available  for  the  low- 
lying  transitions,  allowing  a  validation  of  the  theoretical  methods  for  low 
v,j  states  that  could  then  be  extrapolated  to  high  v,j  states.  This 
extrapolation  could  provide  reliable  data  for  the  prediction  of  parasitic 
lasing  processes  that  can  limit  the  scaling  of  high-power  lasers. 


Av* 
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Abstract 


Analysis  is  presented  of  the  quantum  dynamics  of  the  three-dimensional  F+Hj 
FH+H  reaction  for  total  angular  momentum  J-0.  First,  the  method  (coordinates, 
Hamiltonian,  basis  sets,  close-coupling  method,  and  boundary  conditions)  of 
solving  the  Schrodinger  equation  is  reviewed,  with  emphasis  on  numerical 
construction  of  the  scattering  wavefunction  in  the  region  of  the  collision 
complex.  Then,  four  types  of  analysis  of  the  collision  complex  are  presented: 
(1)  translational  wavefunctions  for  the  dynamically  significant  channels,  (2) 
vibration-rotation  energy  partitioning,  (3)  vibration-rotation  entropies,  (4) 
variation  with  position  along  the  reaction  coordinate  of  the  total  scattering 
wavefunction  density.  Emphasis  is  placed  upon  variations  in  these  quantities 
as  the  system  passes  through  a  quantum  resonance  (near  total  energy  0.36  eV). 
In  part  II  of  this  series,  the  total  scattering  wavefunction  density  and  flux 
are  analyzed  in  the  region  of  the  collision  complex. 
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Three-dimensional  studies  of  this  reaction  are  necessary  to  provide  a 
realistic  basis  for  comparison  with  experiment.  Quantum  studies  using  the 
Born** ,  distorted-wave  Born*^”*\  and  infinite-order-sudden* ^~*^  approxi- 
mations  have  generated  energy-dependent  cross  sections  and  relative 
probabilities  for  forming  HF  products  in  v'j'  states.  Coupled  channel 
studies  (within  the  .^-conserving  approximation)  have  been  presented  by 
Redmon  and  Wyatt.*®”^  Reaction  probability  surfaces  P(E,J)  were  obtained  in 
these  calculations  for  populating  v'a2  and  v'a3  from  H^in  its  ground  state 
(v"0,js0).  These  probability  surfaces  were  summed  over  final  rotational 
states,  and  plotted  as  functions  of  total  angular  momentum  J  and  total 
system  energy  E.  The  v'a2  surface  has  a  maximum  near  0.36  eV  (for  ja0,  hence 
a  classical  impact  parameter  of  b*0  for  ja0).  For  fixed  E  below  this  energy, 
P(E, J)  shows  a  monotonic  decrease  with  increasing  J.  However,  at  energies 
only  slightly  above  0.36  eV,  a  maximum  appears  in  the  probability  cuts  at 
non-zero  J.  This  maximum  appears  at  increasing  values  of  J  as  the  value  of  E 
increases.  The  locus  of  these  probability  maxima  forms  a  three-dimensional 
"resonance  ridge"  on  the  probability  surface.  In  contrast,  the  v»0  v'al 
and  va0  v'*3  reaction  probability  surfaces  show  slow  post-threshold  growth 
with  maximum  reactivity  at  Ja0  for  all  energies  studied  (up  to  0.5  eV).  The 
shapes  of  the  curves  for  higher  J  are  similar,  but  the  threshold  energies  are 
shifted  to  higher  E.  Recently,  Bowman,  Lee  and  Ju  developed  a  method  for 
generating  three-dimensional  probabilities  for  this  reaction  from  collinear 
quantum  probabilities.  These  approximate  probabilities  show  features  in 
good  agreement  with  the  explicit  three-dimensional  calculations,  and  in 
particular,  the  presence  of  the  resonance  ridge  for  the  v'a2  surface.  This 
feature  is  completely  absent  in  results  obtained  within  a  classical 


framework. 
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I.  INTRODUCTION 


During  recent  years,  interest  in  state-to-state  processes  involved  in 
elementary  chemical  reactions  has  continued  to  expand,  and  many  new  and 
informative  results  have  appeared.  The  reaction  F+H^v.j)  FIl(v',j')+H  has 
been  the  reaction  of  greatest  experimental  and  theoretical  effort,  primarily 
because  it  is  a  simple  reaction  with  numerous  dynamical  features  that  can  be 
studied  with  reasonable  accuracy.  On  the  theoretical  side,  quasiclassical 
trajectory  studies  on  a  variety  of  potential  energy  surfaces  have  generated  a 
wealth  of  dynamical  information.  ’  Quantum  collinear  calculations  show  a 
number  of  striking  features  not  predicted  by  the  classical  results. The 
quantum  v“0  v'*2  reaction  probability  curve  exhibits  a  sharp  resonance  peak 

(0.01  eV  half-width)  slightly  above  the  threshold  for  reaction,  while  the  v*0 
v'»3  probability  curve  exhibits  slow  but  monotonic  growth  after  threshold. 
Subsequent  collinear  studies  substantiate  and  extend  these  studies,  *  in 
particular  to  investigations  of  the  effect  of  variations  in  the  topology  of 

Q 

the  potential  energy  surface  on  the  calculated  resonance  features. 

In  a  study  of  the  collinear  reaction,  we  graphically  displayed  the 

scattering  wavefunction,  nodal  patterns,  the  scattering  probability  density, 
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and  the  flux  distribution  as  functions  of  the  reaction  coordinate.  In  that 
study,  we  showed  that  wells  in  the  adiabatic  vibrational  correlation  curves 
and  the  completness  of  the  vibrational  basis  dramatically  influence  the 
resonance  shape.  Recently,  Walker  and  Hayes ^  have  shown  that  the  v“0  v'*2 

collinear  resonance  can  be  destroyed  by  removing  the  entrance  channel 
potential  barrier  (0.05  eV  on  the  Muckermann  surface  V).  These  results 
illustrate  the  importance  of  potential  surface  features  on  these  resonance 
effects,  and  provide  a  measure  of  the  accuracy  desired  in  ab  initio  potential 
energy  surfaces . 


Three-dimensional  studies  of  this  reaction  are  necessary  to  provide  a 
realistic  basis  for  comparison  with  experiment.  Quantum  studies  using  the 
Born**,  distorted-wave  Born*^”*^,  and  inf inite-order-sudden* ^~*7  approxi¬ 
mations  have  generated  energy-dependent  cross  sections  and  relative 
probabilities  for  forming  HF  products  in  v'j'  states.  Coupled  channel 
studies  (within  the  J .-conserving  approximation)  have  been  presented  by 
Redmon  and  Wyatt. *®”^4  geact£on  probability  surfaces  P(E,J)  were  obtained  in 
these  calculations  for  populating  v'm2  and  v'«3  from  t^in  its  ground  state 
(v«0,j*0).  These  probability  surfaces  were  summed  over  final  rotational 
states,  and  plotted  as  functions  of  total  angular  momentum  J  and  total 
system  energy  E.  The  v'm2  surface  has  a  maximum  near  0.36  eV  (for  j*0,  hence 
a  classical  impact  parameter  of  b“0  for  j*0).  For  fixed  E  below  this  energy, 
P(E,J)  shows  a  monotonic  decrease  with  increasing  J.  However,  at  energies 
only  slightly  above  0.36  eV,  a  maximum  appears  in  the  probability  cuts  at 
non-zero  J.  This  maximum  appears  at  increasing  values  of  J  as  the  value  of  E 
increajes.  The  locus  of  these  probability  maxima  forms  a  three-dimensional 
"resonance  ridge"  on  the  probability  surface.  1  In  contrast,  the  v“0  v'*l 
and  va0  v'*3  reaction  probability  surfaces  show  slow  post-threshold  growth 
with  maximum  reactivity  at  J*0  for  all  energies  studied  (up  to  0.5  eV).  The 
shapes  of  the  curves  for  higher  J  are  similar,  but  the  threshold  energies  are 
shifted  to  higher  E.  Recently,  Bowman,  Lee  and  Ju  developed  a  method  for 
generating  three-dimensional  probabilities  for  this  reaction  from  collinear 
quantum  probabilities.  These  approximate  probabilities  show  features  in 
good  agreement  with  the  explicit  three-dimensional  calculations,  and  in 
particular,  the  presence  of  the  resonance  ridge  for  the  v'*2  surface.  This 
feature  is  completely  absent  in  results  obtained  within  a  classical 
framework. 


New  results  from  crossed  molecular  beam  experiments  carried  out  at 
26 

Berkeley,  when  displayed  on  polar  velocity-angle  maps,  support  the 
resonance  interpretation  of  the  quantum  mechanical  calculations.  At  a 
center-of-mass  energy  near  0.08  eV,  the  product  intensity  for  v'=l,2,3 
maximizes  at  0  *  n  (in  the  direction  of  the  incoming  Fluorine  atom). 
However,  when  the  relative  energy  is  increased  to  0.12  eV,  the  v'*2  intensity 
maximizes  near  6  ■  100°,  while  the  v'*l ,3  components  remain  predominately 
backseattered.  This  sideways  peaking  corresponds  to  nonzero  impact  parameter 
collisions,  and  directly  correlates  with  the  quantum  mechanical  prediction 
that  a  resonance  ridge  forms  for  J>0  when  the  collision  energy  exceeds  about 
0.1  eV. 

Differencial  reaction  cross  sections  computed  from  three-dimensional 
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reactive  scattering  in  the  ^-conserving  approximation  agree  very  well  with 

28 

recent  high-resolution  beam  data.  However,  both  the  quantum  mechanical  and 
experimental  results  are  quite  different  than  the  predictions  of  classical 
trajectory  calculations.  Classically,  the  v*=0  v'a2  differential  reaction 
cross  section  remains  backwards  peaked,  even  up  to  0.5  eV  total  energy. 

The  work  presented  in  this  and  the  following  paper  (II)  provides 

the  first  detailed  analysis  of  the  three-dimensional  F+Hj  reaction  in  the 
transition  state  region.  In  addition  to  asymptotic  information  such  as 
reaction  probabilities  and  cross  sections,  we  examine  the  total  scattering 
wavefunction  and  quantities  calculated  with  this  wavefunction,  as  a  function 
of  the  reaction  coordinate.  In  the  present  paper,  by  examining  average 
vibrational  and  rotational  energies,  we  follow  the  partitioning  of  the 
available  energy  as  the  system  proceeds  from  reactants  to  products.  To  gain 
additional  insight  into  the  population  inversion  and  resonance  phenomena,  we 
present  in  II  an  analysis  of  the  scattering  wavefunction  probability  density 
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and  flux  at  three  total  energies:  0.33,  0.36,  and  0.40  eV. 

In  the  following  section  we  discuss  the  coordinates,  Hamiltonian,  and 
the  numerical  solution  of  the  coupled  equations.  In  section  III,  we  discuss 
the  method  used  to  construct  the  three-dimensional  scattering  wavefunction. 
In  section  IV,  the  methods  used  to  analyze  internal  energy  disposal  are 
presented,  and  in  the  final  section  the  insights  gained  into  the  fundamental 
processes  that  dominate  the  dynamics  of  this  important  reaction  are 
summarized. 
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II.  Formulation  of  the  Reactive  Scattering  Problem 
A.  Coordinates 

For  atom-diatom  reactions  at  subdissociative  energies  there  are  three 
arrangement  channels  which  lead  to  separated  A+BC  (channel  a),  AC+B  (product 
channel3  ),  or  AB+C  (product  channel y).  After  removing  the  center-of-mass 
motion,  two  mass-scaled  relative  separation  vectors  describe  the  three-atom 


configuration  in  each  arrangement  channel;  for  example,  in  the  a~channel: 


R  =  <£ 

a  a  c  a  a  b  b  ab 

r  =  c  ( r,  -  r  ) ;  c  =  (u,  /u  ,  ) 
a  a  b  c  a  Mbc  abc 


i 


(D 


Similar  expressions  apply  to  the  other  two  arrangement  channels.  In  these 
coordinates  a  single  reduced  mass 


U  =  M  M.  M  / (M  +M.+M  )  (2) 
abc  abc 

applies  to  all  three  arrangement  channels.  We  now  introduce  a  body-fixed 
frame,  in  each  arrangement  channel  A  ,  with  axes  oriented  relative  to  the 
space-fixed  axes  with  three  Euler  angles  The  BF  frame  is  chosen 
so  that  the  BF  Z-axis  is  parallel  to  for  large  reactant  separations,  but 
changes  smoothly  during  the  course  of  the  reaction  so  that  it  is  parallel  to 
-R  for  the  product  configuration  channel  S  •  This  reorientation  of  the  BF 

p 

axes  is  accomplished  with  a  switching  angle  a(s)  whose  value  depends  upon  the 
extent  of  the  reaction  as  measured  by  the  reaction  coordinate  s.  In  order  to 
specify  the  instantaneous  size  and  shape  of  the  nuclear  triangle  within  the 
BF  frame,  natural  collision  coordinates  (NCC)  are  used.  The 
translational  coordinate  s  varies  from  -°°to  +°°as  the  reaction  progresses. 
In  planes  of  constant  s  perpendicular  to  a  reference  curve  (RC),  usually 
taken  to  be  a  circular  arc  in  the  reaction  zone,  p  and  y  represent 
generalized  vibrational  and  bending  motion.  Asymptotically,  p  and  y 


correspond  to  simple  diatomic  molecule  vibration  and  rotation.  Fig.  1 
illustrates  the  NCC  for  the  a  and  0  channels,  and  shows  a  plane  of  constant 
s. 

Another  useful  set  of  orthogonal  internal  coordinates,  (Z,z,m),  is 

constructed  in  the  following  manner.  By  forcing  the  BF  X-axis  to  be 

perpendicular  to  the  nuclear  plane,  r^  and  have  BF  components  (0,ya»za' 

and  (0,^,za).  Then,  byreplacing  (ya,Ya)with  polar  coordinates  ( m ft), 

where  m  =  r  siny  ,  y  *  — m  coset  ,  and  Y  —  m  sin ct  ,  we  obtain  the  set  (Z,z,m) 

shown  earlier  in  Fig.  1.,.,  The  m-axis  is  perpendicular  to  the  collinear  plane 

(m=0)  and  measures  the  deviation  of  the  nuclear  triangle  from  collinearity . 

Bo  jse  is  a  homonuclear  molecule,  only  two  of  the  three  possible 

arrangement  channels  discussed  above  need  be  considered.  The  a  and  g 

channels  join  at  s*0  (by  definition),  forming  a  surface  on  which  the 
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wavefunction  and  derivative  for  reactant  and  product  channels  must  match. 
The  y  channel  must  also  be  considered  in  completely  specifying  the  boundary 
conditions,  but  not  in  the  integration  of  the  coupled  equations  (for  A+B£ 
systems). 


B.  Hamiltonian  and  Wavefunction  Expansion 

The  NCC  Hamiltonian  operator  can  be  partitioned  into  operators 

representing  pure  translation,  vibration,  and  rotation,  plus  the  remaining 
32 

coupling  terms. 


»  -  «„  *  \ib  *  "rot  *  "coup  <3> 

The  fl  operator  includes  terms  coupling  translation,  vibration,  and 

c^up 

rotation  which  vanish  at  large  s.  Many  of  these  terms  are  small  in  systems 


for  which  the  reaction  passes  through  a  linear  intermediate.  These  terms  are 
neglected  in  the  present  and  previous  NCC  treatments,  and  this  approximation 


□ 


prevents  this  approach  from  consideration  as  an  "exact"  method. 

In  a  similar  manner,  the  electronic  potential  energy  is  decomposed 


V(s,p,y)  ■  V  +  V  ..  +  V.  ,  +  V  (a) 
1  tr  vib  bend  coup 

In  this  representation,  V  is  the  translational  potential  along  the  locus  of 
local  vibrational  minima  (defined  for  straight-line  cuts  perpendicular  to  the 
collinear  RC),  Vv^  is  the  local  vibrational  potential  along  the  cuts,  and 
^ bend  c^e  bending  potential  for  variations  in  Y  away  from  the  local 
vibrational  minima  along  a  circular  arc  of  radius  rQ  in  a  plane  of  constant  s 
(see  Fig.  1) . 

Vibration-rotation  coupling  is  small  for  points  near  the  RP  for  the  present 
surface,  and  has  been  neglected  in  our  work  on  the  F+H^  system.  This 
coupling  could,  of  course,  be  included  without  much  difficulty. 

The  empirical  surface  V  of  Muckermann^  is  used  in  this  and  our 
previous  studies  to  obtain  the  translational  and  vibrational  potentials,  with 
a  simple  functional  form  fit  to  the  M5  bending  potential  for  small  deviations 
from  linearity.  The  bending  potential  is  represented  with  the  function 


Vbend(Y;s)  "  <V0(s)/2)[l  "  cos  2yJ 


where  Vq(s)  is  adjusted  to  represent  the  original  M5  surface  for  small 
deviations  from  linearity. 

The  vibrational  Hamiltonian  is  used  to  construct  locally  adiabatic 
vibrational  eigenvalues  and  basis  functions: 

HvibXv(0js)  *  Ev(s)Xv(p;s)  (6) 


8 


The  vibrational  potential  in  Eq.  4  is  usually  fit  to  a  Morse  potential  so 

that  the  basis  functions  XV(P»*)  are  Morse  oscillators  and  Ey  are  Morse 

energies.  The  lowest  twelve  vibrational  energies  and  the  sums  Ey(s)  +  V£r(s) 

are  shown  in  Fig.  2.  The  wells  in  these  curves  near  8*0  are  produced  by 

expansion  and  contraction  of  the  potential  valley. 

The  rotational  Hamiltonian  is  partitioned  into  terms  representing  a 

symmetric  top,  an  asymmetry  term,  a  Coriolis  term  (coupling  bending  motion  to 
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tumbling  of  the  nuclear  triangle),  and  a  bending  term:  ’ 

*rot  *  Hst  +  Hasym  +  Hcor  +  Hbend  ^ 

Anticipating  the  use  of  locally  adiabatic  rotational  functions,  we  form  a 

symmetric  top  free-rotor  diabatic  basis  by  taking  products  of  normalized 

symmetric  top  functions  and  normalized  associated  Legendre  polynomials, 

J  V 

NMK^X)  Pj(cosy).  These  product  basis  functions  are  eigenfunctions  of 

A  A 

Hst+  H^en<1  (when  0),  for  total  angular  momentum  J  with  projection  K  on 

A  A 

the  BF  Z-axis.  Diagonal izat ion  of  the  representation  of  H  +  H  in  the 

asym  cor 

%  jij 

free-rotor  basis  leads  to  eigenfunctions  (1  that  mix  free-rotor  functions 

j*. 

differing  in  K  but  not  j: 


-JM  V  mJ 

%  ■  L  NHK 


The  free  internal  rotor  function  d>  J a‘  are  Legendre  polynomials  multiplied 
by  terms  dependent  on  £ ,  j,  and  J.^  Asymptotically,  the  are  eigen- 

A 

functions  of  Hrot,  and  j  and  £  are  good  quantum  numbers.  However,  for 

i  a*  Jli 

Vbend  *  ®  t*ie  “t***  representation  of  vj,en<j  in  the  Q  ^  basis  must  be 
diagonalized.  The  equation 


“ro t^j£ ( 9<t>xy  5  8 )  "  *j£ ( 9  ^ j  1 9*yY  * 8  > 


The  adiabatic 


defines  the  hindered  asymmetric  top  functions  and  energies. 
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hindered  asymmetric  top  states  are 


ft 

yv  J 


=  £  <K 


(10) 


JK 

where  R^(y;s)  are  hindered  internal  rotor  functions.  The  j  and  l  values 
no  longer  remain  good  quantum  numbers,  but  are  used  to  label  the  rotor  states 
at  all  values  of  s. 

In  the  present  study,  analysis  is  restricted  to  collisions  with  total 
angular  momentum  J*0.  Thus,  the  orbital  angular  momentum  has  only  the  value 
Jl*j.  Since  M-K-0  for  these  collisions,  the  symmetric  top  functions  are 
independent  of  the  orientation  of  the  three-body  triangle.  Thus;  the  local 

4^  T\i 

asymmetric  top  function  ft  (84>XY»»)  depends  only  on  the  bending  angle  Y,  and 


is  proportional  to  Rjj(Y;s).  It  can  therefore  be  represented  as  a  simple 
expansion  over  Legendre  polynomials 


Rj(Y?s)  =  22  Cj'j^  Pj»(COB  Y) 


(11) 


The  total  scattering  wavefunction  for  input  channel  vij  can  be 
expanded  as: 

<y*(0W  =  (n*/w*r)  22  4'y V  vjf,  (s)  xv'(p;s)  V(yss)  (12) 

v»  J  »«.' 

The  coordinate  dependent  scale  factor  r|^/u)^r  simplifies  the  coupled 
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equations.  The  factor  n  preserves  Hermit ic ityJ ' ,  and  the  factory-5  equals 

the  scaled  atom-molecule  aeparation  for  large  s.  The  expansion  coefficients 

«  t^'e  translational  wavefunctions  (TW),  measure  the  amplitude 


for  occupying  state  v'j' S  '  when  the  reaction  was  initiated  in  state 


C.  Solution  of  the  Coupled  Equations 

Having  defined  the  expansion  of  the  wavefunction,  the  next  step  is  to 

numerically  integrate  a  set  of  coupled  differential  equations  for  the  TV).  By 

the  usual  separation  of  coordinates  into  translational  s  and  internal 

degrees  of  freedom,  we  can  obtain  the  close-coupled  equations  for  the 
J 

f  »./i . o ( a ) .  The  set  of  N  simultaneous  second-order  differential  equations 

V  J  X#  V  J  a* 

is  solved  by  assuming  that,  within  a  small  interval,  the  locally  adiabatic 
basis  sets  are  independent  of  s;  this  eliminates  from  explicit  consideration 
the  first-derivative  term  in  the  equations.  The  actual  distortions  of  the 
basis  sets  with  s  are  included  by  suitable  transformations  at  the  boundaries 
between  adjacent  intervals,  or  "sectors".  Within  one  sector  the  coupled 
equations  are 


al  /  V*  /  / 

,  2  Vj'A'vji,  v"j"invj«, 

as  v"j"£" 

where  the  coupling  matrix  has  elements 


(13) 


DvJ£v'j»*«  =  +  -  (VEv»)/2 

1  *  L  V  (14) 

-  V<3>j<*yln2lxv,>} 

note  that  the  couping  matrix  is  diagonal  in  the  rotational  indices. 
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The  solution  of  Eq.  (13)  via  the  R-matrix  method  ,  followed  by  the 
determination  of  the  scattering  matrix  and  the  TW  can  be  viewed  as  a  three- 
step  procedure.  First,  one  solves  Eq.  (13)  as  a  boundary-value  problem  (R- 
matrix  method).  Second,  asymptotic  boundary  conditions  are  applied  to  define 
the  scattering  matrix  S.  Finally,  the  local  TW  are  constructed  from  the 
correct  asymptotic  conditions  defining  the  appropriate  scattering  states,  by 


11 


repeating  the  propagation  with  appropriate  initial  conditions  using 
propagators  assembled  and  saved  during  the  R-matrix  integration.  This 
process  is  a  new  application  of  the  R-matrix  method  to  reactive  scattering, 
and  details  will  be  more  fully  presented  in  this  and  the  following  sections. 

In  discussing  the  method  of  solution  of  the  coupled  equations  within 
an  arrangement  channel  (tube),  we  will  focus  on  the  even-j  channels  in  the 
reactant  ot  -tube.  The  a -tube  is  partitioned  into  two  regions:  one, 
represented  in  polar  coordinates,  lies  between  s“0  (which  locates  the  a-B 
tube  match  surface  )  and  a  cartesian  (large-s)  region  that  is  separated  from 
the  interior  polar  region  by  another  match  surface,  the  polar-cartesian 
boundary  (PCQ) .  These  two  regions  are  further  subdivided  into  small  sectors 
of  width  h^.  Sector  1  lies  between  Sq  (s-0)  and  Sj,  sector  2  between  s^  and 
s2,  ...  sector  m  between  s^j  and  ao  ,  with  s^  taken  to  be  some  sufficiently 
large  value  of  s.  At  the  center  of  sector  i,  a  unitary  transformation  of 
Outproduces  a  diagonal  energy-dependent  eigenvalue  matrix  A2^(E): 

^(i)g(i)j(i)  _  (15) 

The  coupled  equations  in  the  new  (adiabatic)  representation  can  now  be 
written  as  2N  first-order  equations 


where  the  adiabatic  TW,  £,  are  related  to  the  "diabatic"  TO,  £_,  through  the 
equat ion 

*(i)  -  T(i)  f(i)  (17) 

A  Taylor  series  expansion  of  g,  on  the  LKS  of  sector  i  in  terms  of  &  on  the 
RHS  of  the  sector  permits  us  to  rewrite  Gq.  (16)  as 


The  local  propagators  Pq  are  given  by  analytical  formulas  involving  the 
sector  width  h  and  the  eigenvalues  \  . 

Continuity  of  the  total  wavefunction  and  its  derivative  along  s 
requires  that  at  sector  boundaries  the  adiabatic  TV  meet  the  following 
conditions: 

(jg(i-'l)_  ^(i-1  )g(i-1  ,i)^,(i)(jg(i)  _  (i) 

ar  "  "  =  =  3sl 

The  diabatic  overlap  0ll_  * has  elements  generated  from  overlap  integrals 
involving  the  basis  functions, 


*R 


=  .(1-1  .DJi) 


2.L 


=  F' 


II 


(19) 


i-l.i) 

'J* 


<X 


V  Rj' 


x!iS(i)> 


(20) 


Although  the  adiabatic  TV  are  uncoupled  within  each  sector,  they  mix  at  the 
sector  boundaries. 

Conservation  of  flux  requires  that  the  matrix  product 

T< i  connecting  adjacent  sectors  be  orthogonal;  in  practice  Q 

is  not  orthogonal,  but  it  is  still  diagonally  dominant.  To  insure  the 
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required  orthogonality,  we  use  a  Gram-Schmidt  procedure  at  each  sector 
boundary.  The  combination  of  propagation  across  a  sector  with  matching  at 
the  sector  boundary  can  be  expressed  as  a  matrix  equation  that  integrates  the 
coupled  equations  from  the  RHS  of  sector  i  to  the  RHS  of  sector  i-1 : 


■  Su-1>- 

r 

«f 

■m  m 

R  *- 

,(i-1.i) 


'p(i) 


P 
=  2 


‘  i(1)' 

£{i> 

J 

3s  . 
■»  *■ 

(21) 


In  addition  to  the  sector  matching  requirements,  there  are  three 


additional  matchings,  in  addition  to  the  arrangement  channel  matching,  that 


3 


need  to  be  mentioned.  As  discussed  elsewhere,'"'  at  the  boundary  between  the 
polar  and  cartesian  (asymptotic)  regions,  a  transformation  accounting  for  the 
differences  in  coordinates  must  be  performed.  In  addition,  at  s=0,  the 
matrix  connecting  the  reactant  and  product  channel  adiabatic  TW  is  set 
equal  to  1^ ^  in  order  to  initiate  the  propagation  sequence  at  the  LKS  of 
reactant  sector  1.  (In  our  notation,  the  "zero"  sector  in  the  reactant 
channel  corresponds  to  sector  1  in  the  product  channel).  This  means  that  we 
begin  with  diabatic  TW  at  the  match  surface,  which  facilitates  application  of 
the  matching  conditions.  The  0^’^  matrix  connecting  the  diabatic  TW  at  the 
matching  surface  is  neglected  at  this  step,  but  is  incorporated  into  the  S- 
matrix  equations.  The  third  matching  is  required  to  map  the  scattering 
wavefunction  onto  the  asymptotic  channels,  and  is  incorporated  explicitly  by 
propagation  to  sufficiently  large  s  such  that  all  off-diagonal  elements  of  D 
vanish  and  the  diagonal  elements  become  independent  of  s;  then  ■  f^ra^ . 

Assembly  of  the  local  propagators  and  F  matrices  advances  the 
diabatic  jf  at  s**0  to  the  asymptotic  £  at  8*sm«  This  arrangement  channel 
propagator  is  formed  by  sequentially  multiplying  the  F  and  P  matrices: 


l  E<1>  a<-> 

ds 


Although  this  accumulation  of  local  propagators  is  straightforward,  this 
method  suffers  from  numerical  instabilities  due  to  the  closed  channels  that 
must  be  included  in  the  basis.  The  R-matrix  integration  method  can  be  used 


to  circumvent  these  problems,  because  it  is  inherently  stable  to  the 


inclusion  of  closed  channels. 


V  •« 


The  R-matrix  propagation  equations  may  be  obtained  by  rearrangement 


of  the  propagator  equations'3  into  the  form, 


(i-1) 


1  N11  Kai>ir  J 

=  r(l)  r(1>  -< 

L=3  =4  J  L- 


(i-D- 


The  r^3-^  matrices  contain  matrix  products  of  the  P^3-^  and  and  it  is 

=p  —  q  = 

easy  to  show  that  and  are  symmetric,  and  that  r^  -  r^,  provided 

p(i-  i>i)  £s  orthogonal.  Collecting  the  local  r-matrices  into  an  arrangement 

channel  R-matrix  is  accomplished  not  by  successive  matrix  multiplications,  as 

with  the  propagator  method,  but  by  recursion  formulas  developed  by  Zvijac  and 

Light The  r^3^  are  accumulated  for  each  sector  beginning  with  sector  1  and 
=P 

ending  with  sector  m.  At  sector  m  we  have  obtained  an  arrangement  channel  R- 
matrix  that  relates  the  TW  f  and  their  derivatives  at  s-0  and  s-s  : 


a  ds 


When  these  R-matrices  are  obtained  for  each  arrangement  channel,  they  are 
combined  with  the  appropriate  asymptotic  boundary  conditions  to  allow 
determination  of  the  scattering  matrix.  For  the  initial  channel  v,j,  the 


boundary  conditions  take  the  form: 


Vs)-- k!jr 


[vkv]’ 


•s)  -  S^h^k^s) 
vj  l  vj 


For  the  inelastic  channels,  the  appropriate  form  is 


f  (a)  =  s(-)  s) 

v'j»va;  v*j»  v»j'  x,'  v' j ,8 

while  for  the  reactive  product  channels  the  form  is 
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The  Riccati-Hankel  functions  have  the  asymptotic  forms 


expUOc;:' 

s)  *  ,  lk(i) 

(_expf-[k 


^  exp{-|lr"'s|  -  il-n/2],  %j<0 

while  the  channel  wavenumbers  are  expressed  as 

-  fe  (E  -  E,(±”>  -  -  Vtr(t“>)]  <29> 

The  superscripts  (+)  and  (-)  denote  reactive  and  nonreactive  scattering, 
respectively. 

In  summary,  the  coupled  equations  are  integrated  with  the  R-matrix 
propagation  technique  to  obtain  an  S-matrix.  Then,  a  second  propagation  is 
performed  to  obtain  the  desired  TW  by  using  local  propagators  saved  during 
the  initial  integration.  These  TW  are  used  in  the  analysis  of  the  reacton 
dynamics  discussed  in  the  remaining  sections  of  this  paper  and  in  II. 
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III.  The  Scattering  Wavefunction 

A.  Translational  Wavefunction  Assembly 

In  this  section  we  consider  the  propagation  of  the  TW  and  assembly  of 
the  3-D  wavefunction  at  each  step.  To  initiate  propagation  of  the  TO  in  the 
asymptotic  reactant  or  product  regions,  the  correct  initial  amplitudes  are 
obtained  by  substituting  the  appropriate  S-matrix  elements  into  Eqs.  (25) 
through  (27).  With  these  initial  condition* ,  the  propagation  sequence  in 
Eq.(22)  is  used  to  generate  the  adiabatic  TW  at  each  s,  using  local 
information  saved  at  each  s  during  the  initial  R-matrix  integration.  Eq.(17) 
is  then  used  to  convert  the  adiabatic  &  into  the  diabatic  f.  For  example,  at 
the  mt*1  sector  we  begin  with  the  equation  for  the  even-j  g_  in  the  a  -tube 
(odd-j  construction  proceeds  in  an  analogous  manner): 


R 

rF(m-1,m)p(m) 

F(m-1,m)p(in) 


p(o-1  ,m)p(m)-,  p_(m) 


F(m-1 fm)p(m) 


g 


,  (m) 


(30) 


which  takes  the  adiabatic  TW  and  derivatives  in  sector  m  and  propagates  them 
to  sector  m-1 ,  where  they  are  transformed  into  the  diabatic  representation: 


f(m-l)  _  T(m-1)  g(ra-l)  f , (m-1 )  _  T(m-l)  ,  (m-1)  ^-\) 

Similarly,  to  calculate  g^m”^and  £*  ^  in  sector  (m-2)  we  simply  relabel 
Eq.(30)  by  changing  the  superscripts  m-1  to  m-2  and  m  to  m-1,  substitute  in 
the  and  jg' just  obtained  with  Eq.00),  and  perform  the  matrix 
multiplication.  Successive  application  of  this  procedure  generates  both  the 
adiabatic  and  diabatic  TW  throughout  all  arrangement  tubes.  If  the  numerical 
technique  remains  stable  (remember,  this  is  now  an  initial  value  method), 
1(0)  and  l'(0)  will  satisfy  the  proper  rotor-partitioning  (bifurcation) 


conditions  at  s=0. 


In  the  present  applications  to  F+t^  linear  dependencies  did  arise  in 

the  reactant  tube  in  attempts  to  integrate  large  basis  sets  inward  form  large 

values  of  -s.  The  source  of  the  linear  dependence  of  the  columns  of  £  are 

the  highly  closed  channels,  with  corresponding  small  asymptotic  amplitudes 

and  large  propagator  elements.  Due  to  the  finite  machine  precision,  the 

matrix  multiplications  at  each  step  introduce  local  errors  that  grow  and 

eventually  make  the  value  of  £  meaningless .  Because  of  this,  the  vibrational 

basis  was  limited  to  10  states  in  these  calculations,  and  the  reactant  tube 

propagation  was  started  at  a  modest  value  of  -s.  The  R-matrix  propagation 

does  not  suffer  from  this  problem,  and  for  simply  obtaining  an  S-matrix  the 

number  of  vibratonal  states  and  the  range  of  s  need  not  be  restricted. 

For  the  results  reported  here,  a  60  channel  basis  with  10  vibratonal 

states  and  a  rotor  basis  distributed  12/12/12/8/6/2/2/2/2  was  employed.  The 

rotationally  averaged  results  using  this  basis  agree  reasonably  well  with 
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results  obtained  using  much  larger  basis  sets,  and  they  are  adequately 

converged  for  the  graphical  analysis  of  the  present  work.  In  the  actual 

propagation,  the  even- j  and  odd-j  states  are  propagated  separately. 

Propagation  in  reactants  is  initiated  at  -1.6ag;  at  this  value  of  s  the 

interaction  matrix  is  essentially  diagonal. 

B.  Assigning  Quantum  Numbers  to  the  TW 

One  feature  of  the  TV  that  we  have  not  discussed  is  the  labeling  of 

the  rows  of  f.  and  £  by  quantum  number.  Considering  only  J*0  (£,“j),  this 

means  associating  quantum  numbers  v,j  with  the  elements  of  the  TW.  For  the 

diabetic  TW,  f,  the  assignment  is  straightforward;  one  simply  takes  the 

assignment  used  in  filling  the  interaction  matrix  D.  If,  for  example,  (D) 

=  =  mn 

corresponds  to  channel  (vj),  then  (f)  corresponds  to  f  • ,  since  D  retains 

—  IQ  —  V  J  s 
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Che  same  row  and  column  labeling  at  each  s  in  all  Cubes.  For  Che  adiabatic 
j>,  Che  diagonalization  of  D  in  Che  R-maCrix  propagation  results  in  an 
energetic  ordering  of  the  locally  uncoupled  £.  Thus,  even  asymptotically,  it 
is  not  immediately  true  that  (&)m  =  (f)m>  because  the  diagonal  elements  of  D 
are  not  in  energetic  order.  Therefore,  to  assign  quantum  numbers  to  &,  we 
use  information  contained  in  the  transformation  matrix  at  the  last  sector, 
T^.  This  matrix  contains  only  ones  and  zeros,  where  the  relative  position 
of  the  ones  identifies  the  vj  label. 

C.  TW  Phasing  Techniques 

Our  early  attempts  at  generating  £  produced  inconsistently  phased  TW; 
that  is,  from  sector  to  sector  an  element  (f^n  might  change  sign  many  times 
over  a  single  wavelength.  This  arbitrary  sign  switching  was  due  to  the  phase 
of  the  eigenvector  matrix  in  sector  i  being  independent  of  the  phase  in 

sector  i-1 ;  there  was  no  translational  phase  coherence.  To  overcome  this 

* 

difficulty,  a  phasing  procedure  previously  employed  to  phase  rotor 
35 

functions  was  used  to  phase  the  TW.  Since  the  asymptotic  phase  of  each 
eigenvector  is  arbitrary,  we  initially  multiply  each  vector  by  +1  or  -1  so 
that  the  largest  (in  magnitude)  component  of  the  vector  is  positive.  Then  at 
every  subsequent  step,  we  again  require  the  largest  component  of  the  vector 
to  have  the  same  sign  as  in  the  previous  sector,  multiplying  the  vector  by  -1 
as  necessary  to  achieve  this.  In  phasing  the  rotor  functions  used  in 
assembling  the  TW,  a  similar  procedure  is  employed,  but  account  is  taken  of 
the  asymptotic  phasing  of  Rj  being  not  arbitrary,  but  by  the  fact  that 
Rj(y;s) >Pj(y) ;  thus  the  columns  of  the  rotor  eigenvectors  are  multiplied  by 
-1  as  necessary  to  achieve  this.  This  phasing  is  not  necessary  during  the  R- 
matrix  propagation  since  all  information  on  absolute  phasing  is  lost  in 
forming  the  rotor  overlaps. 
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IV.  Translational  Wavefunction  Analysis  and  Internal  Energy  Disposal 
A.  Diabatic  Translational  Wavefunctions 

This  section  begins  our  analysis  of  the  scattering  wavefunction  by 
first  examining  the  diabatic  TW  and  then  discussing  several 
quantities  that  illustrate  the  partitioning  of  internal  energy  among  the 
channels  during  the  reaction.  Recall  that  the  TW  amplitudes  squared  are 
interpreted  as  local  measures  of  the  probability  density  for  the  individual 
channels.  The  TW  can  reveal  many  detailed  facets  of  the  collision  process 
and  they  serve  as  a  basis  for  understanding  features  of  the  full  3-D 
wavefunction.  Of  the  60  channels  in  these  calculations,  we  will  present  the 
real  part  of  ,f  for  j=0,2,4,  and  6  in  the  three  lowest  vibrational  levels  for 
each  tube.  However,  we  should  mention  that  the  TW  amplitudes  in  the  closed 
vibrational  levels  do  reach  moderate  magnitudes  in  the  reaction  zone, 
although  they  rapidly  decay  as  the  system  moves  out  of  this  region. 

We  first  examine  the  a  -tube  TW,  which  are  shown  in  Fig.  3(a-c)  at 
three  total  energies.  For  purposes  of  illustration,  we  have  scaled  the  a  - 
tube  TW  f^Q  by  0.25  and  f^2  by  0.50.  This  scaling  amplifies  the  contrast 
between  these  two  dominant  channels  and  the  rest  of  the  reactant  tube 
functions.  There  are  several  characteristics  common  to  the  TW  for  all  three 
energies  shown.  With  input  flux  in  f^g  and  with  ISqqI  ■  0.40,  we  expect 
this  channel  to  have  the  largest  amplitude  of  all.  This  input  channel 
couples  strongest  with  f^2  and  consequently  fgj  has  the  next  largest  average 
amplitude.  Within  a  given  vibrational  level  in  the  a  -tube,  the  TW  tend  to 
decrease  in  magnitude  as  j  increases,  and  as  v  increases  from  0  to  2  the 
amplitude  of  a  given  j  level  becomes  smaller.  Since  only  the  j=0  and  2  states 
are  asymptotically  open  in  reactants,  the  wavefunction  will  show  interference 


effects  caused  by  coupling  in  these  two  channels,  the  other  TW  amplitudes 
becoming  vanishingly  small. 


The  next  region  we  will  consider  is  the  match  surface  between  the 
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a  and  p  tubes.  The  rotor  bifurcation  constraints  require  that  at  s=0  the 
a -tube  fvj  split  into  the  $  -tube  f®j+^  and  f^vj»  with  =  f^j ,  and 

Qt  0  Q 

fvj  =  fjj+1  +  fyj.  Comparing  Figs.  3a,  3b,  and  3  c  with  3d,3e,  and  3f 

respectively  confirms  that  these  conditions  are  satisfied.  We  do  not  present 

8  8 

any  of  the  fj\  +  ^  functions  because  they  appear  nearly  identical  to  the  fj\ 

which  are  plotted  in  Fig. 3.  Because  of  the  difference  in  axis  scales  between 
the  a  and  $  tube  plots  and  because  of  the  small  size  of  some  of  the  TW  at 
ss0,  it  is  difficult  to  visually  check  the  accuracy  of  the  amplitude 
partitioning;  numerically  the  agreement  is  to  one  part  in  ten  thousand.  Note 
that  the  phasing  remains  smooth  across  the  s~0  match  surface. 

It  is  in  the  8  -tube  that  we  find  dramatic  evidence  of  the 
vibrational  population  inversion  and  rotational  excitation  of  the  HF  product 
molecules.  At  the  three  energies  studied  here,  all  of  the  channels  shown  in 

Q 

Fig.6(d-f)  are  asymptotically  open,  so  that  at  large  H-HF  separations  the  f^\ 
wavelengths  vary  inversely  with  the  channel  translational  wavenumber> 
Focusing  on  Fig.  3e  for  the  resonance  energy  0.36  eV,  we  can  see  how 
dramatically  the  fgg  amplitude  declines  from  the  very  large  value  it  has  in 
reactants,  so  that  by  s=+1.5  Sq  it  is  negligable.  Similar  comments  hold  for 
the  other  v*0  rotor  states. 


Examination  of  the  v*l  manifold  in  Fig.  3e,we  see  the  beginning  of 

the  extensive  rotational  and  vibrational  excitation  of  the  product  HF 

8 

molecule.  It  is  f^g  that  shows  the  largest  excitation  in  this  figure.  The 
lower  rotors  in  this  manifold  build  amplitude  in  the  region  between  s**+0.5  and 


that  all  of  these  channels  are  energetically  open  in  KF  at  these  total 

go  g 

energies.  An  interesting  feature  of  fj^,  f^  and  fj6  is  that  they  suddenly 
switch  on  at  about  s*0.5.  From  Fig.  3e  we  can  see  the  population  inversion 
process  occuring,  and  we  can  determine  the  FHH  configurations  where  this 
process  begins  to  dominate  the  dynamics  of  the  reaction. 

g 

Comparing  the  functions  channel  by  channel  in  fig.  3(d-f)  we 

observe  that  the  magnitudes  of  the  TV  for  the  0.36  eV  resonant  energy  (again 
interpreted  as  measures  of  the  probability  density)  are  substantially  larger 
than  those  at  0.33  eV,  while  the  j'“4  and  6  functions  have  nearly  double  the 
amplitude  of  the  corresponding  0.40  eV  TV.  Thus,  from  an  analysis  of  the  TV, 
we  have  gained  insight  into  the  characteristics  of  the  vibrational  population 
inversion  process,  and  of  the  amplitude  growth  in  v'**2  at  the  resonance 
energy. 


B.  Vibrational  and  Rotational  Energies  in  the  Collision  Complex 

A  very  interesting  aspect  of  chemical  reactions  concerns  the  disposal 
of  energy  into  the  various  available  channels  as  the  collision  proceeds  from 
reactants  to  products.  Of  particular  interest  are  quantities  such  as  average 
vibrational  and  rotational  energies  and  entropies  as  functions  of  the 
translational  coordinate  s.  One  method  of  calculating  the  average  energies 
is  to  operate  on  the  total  scattering  wavef unction  at  each  s;  however,  this 

g 

direct  approach  produces  equations  that  are  cumbersome.  The  complications  of 
this  approach  arise  because  of  the  curvature  of  the  coordinates.  As  an 
alternative  approach,  we  recast  the  wavefunction  of  Eq.(12)  in  terms  of 
curvature  free  TV  expansion  coefficients  that  eliminate  the  r|^  scale  factor: 

TTj  ■  tE  (32) 

“  r  v'r 
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(33) 


Uv‘j'vj(s)  =  2-J  <Xv''n  6j"j'  fv"j"vj(s) 
v"j" 

The  expansion  of  a  single  IT  mixes  all  of  the  diabatic  f ^  functions  with  the 
same  rotational  index,  each  weighted  by  the  vibrationally  averaged  curvature 
factor.  The  overlap  of  4*  with  itself  becomes 


<We»XOY  ■  iS  (34) 

since  the  3-D  normalization  of  the  vibrational  functions  requires 


v' v 


(35) 


The  U(s)  functions  can  be  utilized  to  define  the  local  normalized 
probabilities  Pyj(s)  *  |  Uvj(s)|  2/£|  Uy j ( s )|  2,  where  the  channels  are  again 
labeled  by  v j .  This  set  of  Pvj(s)  allows  us  to  obtain  average  values  without 
explicitly  operating  on  the  wavefunction.  Thus  the  average  vibrational 
energy  is  obtained  from  the  expression 


<Ev>  =  E  Ev(s)  Pvj(s)  (36) 

VJ 

and  the  average  rovibrational  energy  is  obtained  from 

<Ev  +  Er>  -  E  (vs)  +  wj(s))  Pvj(s)  (37) 

Figure  4  illustrates  these  quantities  for  total  energies 
0.33,  0.36,  and  0.40  eV. 

Looking  first  at  the  average  vibrational  energies  in  Fig.  4a,  we 
notice  that  in  the  asymptotic  reactant  tube  the  average  vibrational  energy  is 
predominately  in  v*0,  as  expected  since  it  is  the  only  open  vibrational 
level.  Approaching  s*0,  the  average  energies  oscillate,  with  the  frequency 
increasing  as  the  total  energy  increases.  As  the  vibrational  energies  plunge 


in  energy  just  before  s*0,  so  do  the  average  vibrational  energies.  However, 
the  average  vibrational  energy  grows  relative  to  the  bottom  of  the  potential 
valley  (the  translational  potential  vCr(s)  has  been  added  to  the  average 
vibrational  and  rovibrational  energies  before  plotting).  Between  s*0  and 
s*+2.0,  the  oscillations  in  the  average  energies  still  reflect  the  mixing  of 
the  upper  vibrational  levels,  but  in  contrast  to  the  reactant  tube,  at  all 
three  energies  the  oscillating  pattern  is  similar;  one  minimum  occurs  very 
close  to  s*0,  another  minimum  occurs  at  approximately  s*+0.4,  and  a  final 
minimum  occurs  near  s**+1.0.  The  amplitude  of  these  oscillations  is  however 
quite  different  for  the  three  levels  because  of  the  increase  in  the  v'*3 
transition  probability  as  the  total  energy  rises  from  0.33  to  0.40  eV.  This 
is  why  at  0.33  eV  the  asymptotic  product  average  vibrational  energy  lies 
close  to  the  v*2  energy  level,  while  at  0.36  eV  more  v'*3  component  mixes  in, 
so  that  by  0.40  eV  the  contribution  from  v'*3  is  sufficiently  large  that 
(  >  is  approximately  halfway  between  the  v“2  and  v-3  levels. 

The  average  rovibrational  energy  plotted  in  Fig.  4b  combines  the 
features  contained  in  Fig.  4 a  with  the  effects  of  local  rotational 
excitation.  These  curves  have  a  similar  interpretation  to  those  in  Fig.  4a, 
with  the  rotational  effects  most  evident  in  the  reactant  tube  (due  to  the 
mixing  of  the  two  open  rotor  states).  An  interesting  feature  of  <E  +  E_>  , 
excitation  into  the  closed  channels  just  prior  to  s*0,  is  indicative  of  the 
quantum  nature  of  this  reaction.  Here,  < Ey  +  ER>  rises  above  the  system 
energy.  The  main  contributor  to  the  tunneling  is  <Ev>.  This  is  not 
surprising  since  quantum  effects  dominate  the  collinear  reaction  as  well. 
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C.  Vibrational  and  Rotational  Entropies  in  the  Collision  Complex 

By  studying  the  local  entropies,  we  can  gain  a  better  understanding 
of  the  distribution  of  probabilities  and  of  the  rovibrational  disorder  during 
the  reaction.  The  total  entropy  at  location  s  (in  units  of  Boltzmanns 
constant,  k)  is  defined  by: 

«,(■)•-£  Pvj(s)  *n£vJ(s))  (38) 

vj 

Two  other  entropies,  the  rotationally  averaged  vibrational  entropy,  R, 

and  the  vibrationally  averaged  rotational  entropy,  < y  require  averaging 
the  Pvj  functions  in  the  following  manner: 

<SV>R  =  PVj)  Pvj) 

<vv  =  *(£%,).  . 

j  V  /  V 

Figure  5  illustrates  these  three  entropies  at  the  energies 
0.33,  0.36,  and  0.40  eV.  Since  each  type  of  entropy  for  these  three  energies 
shows  similar  behavior,  we  will  not  consider  their  energy  dependence,  but 
instead  will  emphasize  the  general,  features  of  each  entropy. 

Looking  first  at  (Sy>R,  we  see  that  in  asymptotic  reactants  there 
is  zero  entropy  because  only  v*0  is  open.  But  as  s~0  is  approached,  this 
entropy  rapidly  grows  and  then  falls  off  to  a  nonzero  value  for  s>1.5.  This 
shows  both  the  short-lived  excitation  in  the  transition  state  due  to  the 
participation  in  the  scattering  of  asymptotically  closed  vibrational  levels, 
and  the  uncertainty  of  the  final  vibrational  state  of  products.  The  nonzero 
value  at  large  product  separation  simply  reflects  the  partitioning  of  amplitude 
between  v'*2  and  v'»3.  Considering  next  the  <SR>y  plot,  we  notice  that  in 
reactants  there  is  a  small  (compared  to  products)  rotational  entropy  due  to 
mixing  between  j'«0  and  j'*2.  Note  that  just  before  s>0  the  entropy  curves 
nearly  reach  zero.  Referring  back  to  Fig. 3  ,  it  can  be  seen  that  almost  all 


25 


the  flux  flows  through  the  lowest  (j'*0)  adiabatic  bending  channel  in  this 


region.  However,  due  to  the  large  number  of  available  open  rotational 
channels  in  products,  y  becomes  quite  large  for  s>1.0. 

The  curves  contain  features  of  both  <Sy>  R  and  <SR>  v:  nonzero 
entropy  in  reactants  and  products  with  the  entropy  in  products  greater  than 
in  reactants,  and  a  peak  near  s*0. 

This  concludes  our  analysis  of  the  translational  wavef unctions  and 
their  associated  average  energies  and  entropies.  We  have  focused  attention 
on  the  flow  of  amplitude  through  the  individual  channels  as  functons  of  both 
the  system  energy  and  reaction  coordinate. 


D.  Integrated  Density  vs  Reaction  Coordinate 

We  acquire  an  additional  perspective  on  the  collision  process  by 
compressing  information  from  the  wavefunction  probability  density. 
Integration  of  the  total  density  within  planes  of  constant  s  then  allows 
construction  of  an  interesting  plot  of  the  integrated  density  as  a  function 
of  the  translational  coordinate  s.  The  integrated  density  is  defined  by  the 


equat ion 


••// 


¥  (n,m;s)  'P(n,m;s)  dn  dm 


The  factor  b)  appears  in  Eq.  (40)  to  prevent  the  density  from  decaying  asymp- 
iotically  as  1/R  .  Fig. 6,  illustrates  the  s-dependence  of  the  integrated 
density  for  the  energies  0.33,  0.36  and  0.40  eV.  The  oscillations  in  the 
densities  are  caused  by  interference  between  all  locally  open  channels.  The 
frequency  of  the  oscillations  can  be  correlated  to  some  degree  with  the  fact 
that  available  translational  energy  increases  with  increasing  total  energy. 
The  decrease  in  magnitude  of  the  densities  past  about  s*-0 . 1  is  a  consequence 
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of  the  considerable  amount  of  nonreactive  flux  in  this  reaction.  The 
resonant  curve  (the  0.36  eV  density)  reaches  a  maximum  before  s*0,  and  then 
decreases  as  the  system  moves  into  the  product  region.  Thus  we  might  call 
this  v'm 2  resonance  an  "entrance  channel  resonance”.  The  vibrational 
distributions  have  mostly  settled  down  by  8*0.5,  and  only  the  rotational 
populations  show  much  change  after  this  region  of  s. 


V.  Conclusions 


The  quantum  dynamics  of  the  three-dimensional  F+H2  reaction  has  been 
analyzed  for  three  total  energies,  and  for  total  angular  momentum  J*0.  The 
energy  0.36  eV  is  of  special  significance,  since  this  is  close  to  the  J*0 
resonance  energy  for  the  process  F+H2(v-0)-» H+HF(v'-2) .  The  computational 
methods  used  to  obtain  the  scattering  wavefunction  required  for  the  dynamical 
analysis  in  the  region  of  the  collision  complex  have  been  presented.  This 
has  been  prefaced  by  a  brief  review  of  the  coordinates,  Hamiltonian,  basis 
sets,  and  boundary  conditions  used  in  the  NCC  approach  to  the  reaction 
dynamics . 

In  the  dynamical  analysis,  emphasis  was  placed  on  graphical  methods 
for  illustrating  the  variation  with  s  of  the  translational  wavef unctions, 
vibrotor  energies,  vibrational,  rotational  and  total  entropies,  and  the  total 
density  in  planes  perpendicular  to  the  reaction  coordinate.  We  will  briefly 
summarize  some  of  the  results  from  this  analysis.  In  Fig.  6  the 
translational  wavef unctions  fyj(s)  are  shown  for  the  j“0-6  rotational  states 
associated  with  the  v«0-2  vibrational  manifolds.  In  the  F+H2  reactant  tube, 
the  translational  wavefunction  f02  has  an  amplitude  which  is  at  most  one  half 
of  (qq.  Hear  s*0,  there  is  a  small  amount  of  excitation  into  the  lower 
rotational  levels  of  the  vs0  and  1  manifolds.  In  the  H+HF  product 

g 

arrangement,  f^  decays  almost  completely  by  s*+1.0  Sq,  while  amplitude 
builds  up  in  the  j*4-6  rotational  levels  of  the  v*0  and  1  manifolds. 
However,  the  most  dramatic  effect  in  the  product  region  is  the  relatively 
rapid  buildup  between  s*0.5  and  8*1.5  of  amplitude  in  the  j*4-6  rotational 
levels  of  the  v*2  manifold.  This  rotational  excitation  occurs  on  the  final 
part  of  the  turn  into  the  product  valley  (see  Fig.  4). 
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The  oscillations  seen  in  Fig.  7  of  the  average  vibrational  energy 
<Ey>  and  the  average  rovibrational  energy  <  E^  +  ER>  as  s  goes  from 
reactants  to  products  are  associated  with  interference  effects  among  the 
vibrational  channels.  Of  more  interest  than  these  average  energies  are  the 
average  and  total  entropies  shown  in  Fig.  8.  In  reactants,  the  rotationally 
averaged  vibrational  entropy  is  zero  (only  v**0  is  populated),  but  the 
vibrationally  averaged  rotational  entropy  is  nonzero  (j*0-2  are  populated). 
Near  8*0,  the  rotational  entropy  is  lower  than  in  reactants  or  products 
because  of  the  strong  propensity  for  the  lowest  bending  state  of  the  complex. 
In  products,  the  vibrational  entropy  is  about  half  the  peak  value  in  the 
collision  complex,  because  flux  is  leaking  out  of  the  temporarily  excited 
virtual  states  (primarily)  into  the  v*2  manifold.  Also,  the  rotational 
entropy  in  products  is  large  (about  twice  the  vibrational  entropy),  due  to 
the  relatively  diffuse  rotational  excitation  through  many  rotational  channels 
(but  primarily  j*4-6). 

The  variation  with  s  of  the  total  wavef unction  density,  shown  in 
Fig.  9,  measures  the  net  tendency  of  probability  to  accumulate  at  position  s. 
This  quantity,  I(s),  is  computed  as  an  integral  over  a  plane  of  constant  s  of 
the  local  probability  density.  The  most  interesting  features  of  this 
quantity  are  found  to  be:  the  oscillations  in  reactants  due  to  interference 
between  incoming  and  outgoing  waves,  and  the  factor  of  10-100  drop  in  I(s)  on 
turning  the  corner  into  the  product  valley.  However,  when  comparing  plots  of 
I(s)  at  the  three  energies  considered  in  this  paper,  the  region  between 
s*-0.5  &q  (near  the  entrance  channel  barrier)  and  s*0  is  especially 
significant.  On  resonance  (0.36  eV),  the  density  in  this  regon  is  about  a 
factor  of  two  or  more  larger  than  for  the  other  two  energies,  which  lie  above 
and  below  the  resonance  energy.  The  resonance  is  associated  with  buildup  of 


amplitude  between  the  entrance  channel  barrier  and  the  turn  into  products; 
thus  it  is  termed  an  "entrance  channel  resonance". 
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In  part  II  of  this  series  ,  density  and  flux  maps  obtained  from  the 
total  scattering  wavefunction  will  be  shown  in  planes  perpendicular  to  the 
reaction  coordinate  and  in  planes  for  bent  geometries  (m>0)  parallel  to  the 
collinear  plane.  Emphasis  will  be  placed  on  the  behavior  of  these  quantities 
as  the  system  goes  through  resonance. 
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FIGURE  CAPTIONS 


1.  Natural  collision  coordinates  for  3-D  reactive  scattering.  Z^  and  z 
are  the  components  of  the  mass  scaled  F-Hj  and  HH  relative  separations 
on  the  body-fixed  Z-axis.  Z^  and  (not  shown)  are  the  components  of 
the  corresponding  product  separations  on  the  BF  Z-axis.  Nonlinear 
geometries  are  measured  with  m;  the  ms0  plane  is  the  collinear  plane. 
The  reaction  coordinate  s  is  measured  along  the  reference  curve  RC.  In 
the  reactant  and  product  Cartesian  regions,  planes  of  constant  s  lie 
perpendicular  to  Z r  and  Z^,  respectively,  while  in  the  polar  region 
constant  s  planes  intersect  the  collinear  plane  on  lines  connecting  the 
turning  center  TC  and  the  reactant  or  product  Z-axes.  The  floating 
origin^®  FO  for  r  and  y  varies  with  s  and  is  introduced  to  simplify 
terms  in  the  rotational  Hamiltonan.  The  TC  is  located  at 
(7.000,  2.579),  and  the  radius  of  the  circular  reference  curve  is  1.563. 
The  cartesian-polar  boundaries  are  at  s*-1.904  and  s*+1.740  in  reactants 
and  products,  respectively.  At  large  to  moderate  separations,  y  is  the 
geometric  angle  between  the  molecule  and  atom-molecule  vectors. 

2.  Adiabatic  energy  correlation  diagrams  (energies  in  eV) .  (a)  vibrational 

energies  ;  (b)  vibrational  energies  plus  the  translational  potential 
vtr  »  (c)  hindered  asymmetric  top  energies  Wj  for  J“0;  (d)  the  sum  of 
the  translational  potential,  lowest  vibrational  energy,  and  top 
energies. 
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3.  (a)-(f)  The  real  part  of  the  diabatic  translational  wavefunctions  at  the 
total  energies  0.33,  0.36,  and  0.40  eV.  Only  four  vibrotor  TO  in  each 
of  the  v=0,l  and  2  manifolds  for  the  60  channel  basis  are  shown. 
(a)-(c)  are  reactant  tube  TW;  (d)-(f)  are  product  tube  TO. 

4.  (a)  The  average  vibrational  energy  and  (b)  the  average  rovibrational 
energy.  The  average  energy  curves  are  shown  with  dots,  and  the 
background  curves  are  the  vibrational  energy  levels. 

5.  Entropies  as  a  function  of  s  for  the  three  system  energies  considered 
(units  are  k,  Boltzmann's  constant). 

6.  A  semi-log  plot  of  the  integrated  s-plane  wavef unction  densities  for  the 
total  energies  0.33,  0.36,  and  0.40  eV. 


Figure  3c 
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ABSTRACT 

Velocity-scattering  angle  intensity  maps  for  the  F+^Cv-O.j-O)  -»■ 

FH(v'*2, j *)+H  reaction  are  predicted  from  quantum  mechanical  J  -conserving 

z 

calculations.  The  shift  in  the  angular  distribution  from  backscattering 
at  1.8  kcal/mole  to  sideways  scattering  (intensity  peak  at  100°)  at  3.0 
kcal/mole  is  in  quantitative  agreement  with  recent  high  resolution  crossed 
molecular  beam  experiments. 
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Quantum  mechanical  resonance  structure  in  the  collinear^"  F-H^v^O) 

ra(v'*2)+H  reaction  was  predicted  in  1973  and  in  the  three-dimensional 
2 

reaction  in  1979.  In  order  to  illustrate  their  results,  Redmon  and  Wyatt 
plotted  the  reaction  probability  surface  for  v*0  -*■  v'“2  reaction.  This  is 
a  surface  of  total  probability  for  the  reaction  (v*0,j*0)  -*■  (v'*2,Ej'), 
plotted  vs.  total  energy  E  and  total  angular  momentum  J  (since  j*0  in  reac¬ 
tants,  J“il,  where  l  is  the  atom-molecule  relative  orbital  angular  momentum; 
l  is  related  to  the  Impact  parameter,  yvb*»$,h) .  This  reaction  probability 
surface  is  shown  In  Figure  1;  as  E  is  Increased  at  each  value  of  J,  the 
reaction  probability  suddenly  increases,  reaches  a  peak  at  the  resonance 
energy  Er(J),  and  then  declines  onto  a  shoulder  at  higher  energy.  The 
"resonance  ridge"  of  high  reactivity  moves  to  higher  E  as  J  Increases. 

(The  J  dependence  of  Ef  is  remarkably  well  predicted  by  a  new  semiclasslcal 
model  based  upon  resonant  periodic  orbits  ).  Up  to  0.35  eV,  the  reaction 
probability  declines  monotonically  with  J.  However,  for  E<0.37  eV,  the 
peak  reactivity  is  reached  for  a  nonzero  value  of  J  (e.g.,  at  0.40  eV, 
the  peak  reactivity  occurs  at  about  J»10) .  By  way  of  contrast,  the  0-*-3 
reaction  probability  surface  decreases  monotonically  with  increasing  J, 
for  E  £  0.50  eV.  The  purpose  of  this  Letter  is  to  show  that  these  quantum 
mechanical  results  also  predict  that  the  0-*2  differential  reaction  cross 
section  undergoes  a  dramatic  change  when  the  energy  is  Increased  from  0.35 
to  0.40  eV. 

In  their  1979  crossed  molecular  beam  studies  of  the  F-H^  reaction, 

4 

Sparks  et  al.  found  results  which  are  qualitatively  consistent  with  the 
quantum  mechanical  reaction  probability  surfaces.  When  the  product  in- 


tensity  was  plotted  in  terms  of  velocity-scattering  angle  intensity  maps, 
then  at  a  mean  relative  collision  energy  of  2.0  kcal/mole,  the  product 
molecules  for  v'*l,  2,  and  3  were  all  scattered  predominately  in  the 
8  »  180°  backward  direction  (i.e,  back  in  the  direction  from  which  F 
approached  the  Hj  molecule.)  However,  when  the  mean  collision  energy 
was  Increased  slightly  to  3.2  kcal/mole,  the  peak  in  the  v'=*2  product 
distribution  moved  to  a  smaller  angle  —  extensive  sideways  scattering  occurred. 
By  way  of  contrast,  the  v'*l  and  3  distributions  remained  primarily  back- 
scattered  at  the  higher  energy.  On  the  basis  of  the  reaction  probability 
surfaces,  and  the  new  quantum  results  reported  here,  we  Interpret  the 
shift  to  sideways  scattering  in  the  HF(v’“2)  product  distribution  as  a 
direct  manifestation  of  the  quantum  resonance. 

Before  presenting  results  for  the  F-H^  reaction,  we  will  briefly 
review  the  relevant  theory^.  In  an  atom-diatom  collision  at  total  angular 
momentum  J,  if  we  were  following  the  close-coupling  (CC)  procedure,  the 
total  scattering  wavefunction  would  be  expanded  in  products  of  local  vibra¬ 
tional  and  rotational  basis  functions.  In  our  approach,  the  rotational 
basis  functions  are  hindered  asymmetric  top  states^,  which  are  labeled  by 
the  asymptotic  quantum  numbers  j  and  it.  Here  j  is  the  rotational  quantum 
number  of  the  diatom  and  it  denotes  the  atom-molecule  relative  (orbital) 
angular  momentum  quantum  number.  The  triangle  inequality  limits  the  range 
of  allowed  values  for  it:  <  it  <  J+j  (parity  selection  rules  further 

limit  the  number  of  it's  in  this  Interval).  Except  for  very  low  energies, 
the  number  of  vibration-rotation  channels  (vjit)  in  a  dose-coupling  cal¬ 
culation  at  energy  E  and  total  angular  momentum  J  is  so  huge  that  we  are 


prevented  from  obtaining  "exact"  solutions  of  the  CC  equations.  For  this 

reason,  we  are  forced  to  use  a  channel  reduction  scheme  to  cut  the  total 

number  of  channels  down  to  a  workable  level.  In  our  version  of  the  J  -con- 

z 

g 

serving  approximation  ,  which  was  used  in  all  calculations  for  J  >  0,  we  use 
only  one  asymmetric  top  function  in  the  wavefunctlon  expansion  for  each 
value  of  J  and  j.  This  special  value  of  H  (denoted  L)  is  selected  auto- 

g 

statically  such  that  these  asymmetric  top  states  produce  a  nunHimim  concen¬ 
tration  of  amplitude  for  near-collinear  geometries  in  the  collision  com¬ 
plex.  The  procedure  is  called  ".^-conserving"  because  the  asymmetric 
top  states  which  are  used  always  concentrate  amplitude  in  the  FHH  K-0 
bending  states  (K  is  the  angular  momentum  quantum  number  for  twist  about 
the  z-axis  of  the  complex) .  Solution  of  the  reduced  set  of  CC  equations 
then  leads  to  S -matrix  elements  (amplitudes  for  outgoing  waves)  which  are 


labeled  S  .  _  .  ,  or  since  we  are  interested  here  in  reactions  initiated 

VoV^1 

from  the  ground  state  of  HjCv-O.j-O) ,  we  will  use  a  shortened  notation  for 
the  S-matrix  elements:  (where  L  is  determined  uniquely  by  J  and  j). 

Once  the  S-matrix  elements  have  been  computed,  we  can  use  them  to 
construct  both  differential  and  total  reaction  cross  sections.  The 
differential  reaction  cross  section  is  computed  from  the  hellcity  amplitude: 

hvj(6;E)  -  (^’ko)'1Jo(2j+l)ls  *j(9)SvjL(B)» 


2  2 

where  k  is  the  reactant  channel  translational  wavenumber  (Efc  _  ■  h  k  /2 u  . 

o  crans  o  r 

Vr  ■  mA(mj-hn^)/(m^+mB+m(,)),  and  where  Pj(9)  is  a  normalized  Legendre  polynomial. 
The  angle  0  is  defined  so  that  0-180°  corresponds  to  "backscattering"  of  FH 
in  the  direction  from  which  F  approached  the  Hj  molecule.  The  differential 


reaction  cross  section  is  then 


l 
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avj(0;E)  -  |hvj(0;E)|2.  (2) 

Finally,  the  total  reaction  cross  section  Is  given  by 
o  (l)  -  /0”  .In0d0  o  (6,E) 

--?-j„<2J+1)|SiJ2-  <3> 

o 

The  S -matrix  elements  are  complex  numbers  whose  magnitude  and  phase  enters 

the  computation  of  (0;E) ,  while  only  magnitudes  enter  the  total  cross 

section.  Equation  (3)  Is  the  same  one  that  we  previously  used  to  obtain 

2 

total  cross  sections  from  the  S -matrix  . 

9 

The  J  -conserving  calculations  were  performed  at  two  total  energies 
z 

(measured  from  the  bottom  of  the  F+Hj  entrance  valley;  Etrans  ■  E-0.27  eV)  , 

0.35  and  0.40  eV,  which  correspond  to  relative  collision  energies  of  1.8  and 

3.0  kcal/mole,  respectively.  At  each  energy,  J  -conserving  calculations  were 

z 

performed  for  each  value  of  J  between  0  and  18.  The  basis  sets  always  em¬ 
ployed  10  vibrational  manifolds,  with  the  following  rotor  distributions: 
16/14/12/12/6/2/2/2/2/2,  0.35  eV; 

18/18/18/16/14/12/6/4/2/2,  0.40  eV. 

This  notation  means  that  a  total  of  16  hindered  rotor  states  were  used 
in  the  v*0  manifold  at  0.35  eV,  etc.  In  spite  of  the  fact  that  a  larger 
basis  was  used  at  the  higher  energy  (110  channels,  vs.  70  at  the  lower 
energy),  the  results  at  0.40  eV  are  less  converged  than  at  0.35  eV.  Parti¬ 
cularly  for  J  >  10,  the  S-matrices  are  not  fully  converged  with  respect 
to  shifting  rotor  basis  functions  from  one  manifold  to  another.  In  spite 


6 


of  this,  we  believe  that  the  results  presented  here  have  at  least  "semi- 
quantltatlve"  validity. 

It  was  mentioned  earlier  that  the  helicity  amplitude  (and  hence  the 
differential  cross  section)  depends  upon  both  mangitudes  and  phases  of 
Individual  S -matrix  elements.  Lack  of  convergence  in  the  S -matrix  elements 
(particularly  the  phases)  can  lead  to  artificial  wiggles  In  the  computed 
differential  cross  sections.  In  order  to  avoid  artificial  oscillations  in 
the  cross  sections,  we  carefully  examined  the  "smoothness"  of  both  the 


magnitudes  and  phases  of  all  S -matrix  elements,  at  both  energies.  Plots 
were  made  of  vs •  J  at  each  value  of  j  for  v-2 .  In  addition,  if 

denotes  the  phase  of  an  S -matrix  element,  plots  were  made  of  4^^  and 
-  4>J.T )  vs.  J  (again,  for  each  j  at  each  E) .  In  some  cases  (usually 

VJ  L  VJ  L 

higher  J  values),  magnitudes  and  phases  were  adjusted  "by  eye"  to  produce 
smooth  looking  plots.  In  spite  of  the  somewhat  subjective  character  in 
this  procedure,  we  felt  that  it  compensated  to  some  extent  for  convergence 


errors  in  the  S -matrices. 


The  differential  cross  sections  for  the  F+fljCv^O.j^O)  FH(v’«2, j ')+H 
reaction  at  0.35  and  0.40  eV  are  shown  in  Figures  2  and  3,  respectively. 

The  results  are  shown  as  intensities  I(V,6)  on  velocity-angle  maps,  where 
"velocity"  refers  to  the  recoil  speed  of  FH  and  H  in  the  v'j'  final  channel, 
and  where  "angle"  refers  to  the  reactive  scattering  angle  0  in  Eqs.  (1)  and 
(2).  A  contour  map,  with  cross-hatched  maxima,  is  shown  in  the  left  panel 
of  each  figure,  and  a  more  dramatic  perspective  plot  (viewed  from  the 
"forward"  direction)  is  shown  in  the  right  panel.  Of  course,  only  specific 


final  speeds  (radii)  corresponding  to  v,j’>(2,0),  (2,1),  ... 


are  possible. 


but  the  contour  and  perspective  routines  automatically  fill  in  other  values 
to  make  continuous  functions  for  plotting.  Looking  first  at  Figure  2,  the 
intensity  is  clearly  peaked  at  9  ■  180s  (backscattering) ,  and  the  dominant 
rotor  distribution  is  for  j '*4-5 .  There  is  a  smooth,  monotonic  decline  in 
I(V,9)  as  @  moves  to  the  forward  direction;  practically  no  intensity  is 
found  in  the  forward  hemisphere.  Turning  now  to  Figure  3,  we  note  the 
dramatic  qualitative  change  In  the  intensity.  The  intensity  has  now  shifted 
so  that  it  now  peaks  at  9  ■  100°,  a  full  80s  off  the  backward  direction,  with 
I(V,  100°)/1(V,  180°)  ~  2.  In  addition,  the  rotor  distribution  peaks  at  j'“3; 

the  cross-hatched  peaks  lie  at  slightly  smaller  radii  (speeds)  than  the  outer 
limiting  circle  (corresponding  to  v'*2,  j ’-0) .  (The  rotational  distribution 
peaks  at  higher  j'  values  for  calculations  with  other  basis  sets;  this  point 
is  under  investigation).  Also  apparent  in  Figure  3  is  the  small  backward 
peak  at  0  =  180°;  unlike  the  major  peak  near  0  *  100°,  the  height  of  this 
peak  is  particularly  sensitive  to  phases  of  the  S-matrix  elements.  At  this 
point,  we  cannot  say  definitively  that  there  is  really  a  small  ripple  in 
I(V,0)  In  the  backward  direction.  We  are  much  more  confident  in  claiming 
that  the  large  peak  near  100°  is  a  reflection  of  the  dynamics  and  the  under¬ 
lying  potential  surface. 

The  sideways  shift  in  the  v'»2  quantum  mechanical  product  intensity 
map  is  In  agreement  with  intensity  maps  from  crossed  molecular  beam  experi¬ 
ments.  In  recent  high  resolution  F+Hj  experiments^  (which  came  several  years 
after  the  1979  results  that  were  reported  in  ref.  4)  in  which  the  most  probable 
collision  energy  was  3.1  kcal/mole,  with  a  FWHM  of  only  0.04  kcal/mole,  the 
product  Intensity  for  v'*2  peaked  at  about  100°,  with  a  most  probable  rota¬ 
tional  quantum  number  of  about  j’-7.  In  addition,  the  v'»2  intensity  near 
100°  is  about  a  factor  of  1.7  higher  than  in  the  backward  direction.  Al- 


8 


h 


h 


though  our  current  calculations  underestimate  the  extent  of  rotational  exci¬ 
tation  in  the  vf*2  manifold,  the  shift  of  the  intensity  peak  as  the  colli¬ 
sion  energy  is  increased  from  2  to  3  kcal/mole  is  in  quantitative  agreement 
with  the  high  resolution  beam  results! 

Comparisons  of  these  results  with  other  theoretical  predictions  are  of 
Interest.  First,  quasiclass leal  trajectory  calculations  on  the  Muckerman  V 
potential  surface,  even  at  a  collision  energy  of  5  kcal/mole,  fall  to  predict 
sideways  scattering  of  the  HF(v**2)  molecules^.  Classical  mechanics  com-  1 

pletely  misses  the  sideways  shift  in  the  v’*2  Intensity  that  occurs  between 

12 

0.35  and  0.40  eV!  Turning  to  quantum  calculations,  IOS  calculations  have 

been  performed  at  selected  energies  between  0.31  and  0.50  eV  on  the  Muckerman  ^ 

V  potential  surface.  The  v-2  differential  reaction  cross  section  at  0.36  eV 

is  backpeaked,  but  shows  a  slight  shoulder  between  80*  and  120°.  When  the  energy 

is  Increased  to  0.423  eV,  the  differential  reaction  cross  section  has  about  the  same 

values  at  85°  and  180°,  but  is  10%  lower  at  140°,  and  declines  rapidly  for  9  £  70°.  1 

At  0.50  eV,  there  is  a  slight  sideways  peak  at  9  ■  40°;  the  height  at  40°  is 

40%  higher  than  at  180*.  There  is  clearly  a  trend  toward  sideways  scattering 

at  higher  energies  in  the  IOS  results,  but  the  extent  of  the  angular  shift 

between  0.35  and  0.40  eV  that  appears  in  both  the  J  -conserving  calculations, 

z 

and  in  the  experimental  results,  is  not  quantitatively  predicted.  This  is 

probably  due  to  an  important  difference  between  the  J  -conserving  and  the 

IOS  reaction  probabilities;  in  the  ^-conserving  case,  for  E  2  0.37  eV,  the 

peak  reaction  probability  in  v'"2  occurs  for  J  >  0  (resonance  ridge  in  Figure  1), 

while  the  IOS  calculations  predict  flat  reaction  probabilities  between  J*0  and 

Jnax  “  10  at  0.423  eV) .  Another  approximate  quantum  procedure  for  gen- 
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eratlng  differential  cross  sections  has  recently  been  proposed  .  S-matrix 
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elements  for  v+v'  colllnear  reactiona  were  used  (at  shifted  energies  to 

account  for  both  centrifugal  effects  and  the  bending  energy  of  the  complex) 

In  Eq.  (2)  to  generate  "rotationally  summed"  cross  sections.  In  going  from 

Etrans  “  up  to  2,9  kcal/mole,  shifts  from  backpeaked  to  sidepeaked, 

with  a  maximum  near  165°.  In  addition,  the  ratio  is 

about  2.5  at  the  higher  energy.  In  this  model,  the  0  +  2  reaction  probability, 

when  plotted  vs.  J,  peaks  at  J«*12  at  2.9  kcal/mole,  but  monotonically  declines 

from  its  value  of  J-0  at  the  lower  energy.  These  trends  are  in  qualitative 

agreement  with  our  J-conserving  results.  Finally,  DWBA  differential  reac- 

z 
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tion  cross  sections  have  been  computed  at  several  total  energies  between 

0.31  and  0.57  eV.  Comparing  results  at  0.37  eV  and  0.57  eV,  the  peak  in 

a0+2  moves  f°rwar<*  from  180°  to  110°.  Again,  we  see  a  trend  to  sidepeaking, 

but  neither  the  onset  nor  the  extent  of  sideways  scattering  is  as  abrupt  as 

in  both  the  J-conserving  results  and  the  experimental  beam  results, 
z 

In  summary,  from  quantum  mechanical  calculations  within  the  J-conserving 
approximation,  product  Intensity  maps  have  been  generated  for  the  F+J^Cv-O.j-O)  + 
FH(v'-2,j ')+H  three-dimensional  reaction.  The  extent  of  the  sideways  shift  in 
the  angular  distribution  is  in  quantitative  agreement  with  recent  high  reso¬ 
lution  crossed  molecular  beam  results.  We  believe  that  this  sideways  shift 
is  a  direct  manifestation  of  a  quantum  mechanical  resonance  in  the  FHH  colli¬ 
sion  complex.  The  nature  of  this  resonance  (including  scattering  wavefunctlon 
density  and  flux  maps  in  the  transition  state  region,  entropy  and  vibrotatlonal 
energy  analysis  along  the  reaction  coordinate,  Argand  diagrams  showing  the 
energy  dependence  of  selected  S-matrix  elements,  phase  shift  analysis,  time- 
delays  computed  from  the  energy  derivatives  of  the  S-matrix  phases,  and  a 
kinematic  model  for  the  sideways  shift  in  the  angular  distribution)  is  explored 
in  more  detail  elsewhere^. 


1.  Quantum  mechanical  reaction  probability  surface  for  the  F+H(v«0, 
j«0)  ■+  FH(v'“2,Zj ')+H  reaction  (adapted  from  ref.  2).  J  is  the  same  as 
l,  the  atom-molecule  orbital  angular  momentum  quantum  number.  E„_  is  the 
total  energy,  measured  from  the  bottom  of  the  entrance  valley. 


2.  Product  intensity  for  v'*2  in  velocity-scattering  angle  variables.. 
The  energy  is  0.35  eV  (the  total  v«0  -*■  v'*2  reaction  cross  section  is  1.2a  ;. 
The  left  panel  shows  a  contour  map  of  the  intensity,  with  cross  hatched  ° 
maxima.  The  right  panel  shows  a  perspective  view  of  the  intensity,  viewed 
from  45°  off  the  forward  direction. 


3.  Product  intensity  for  v'*2  in  velocity-scattering  angle  coordinates. 
The  energy  is  0.40  eV.  At  this  energy,  the  total  v-0  -*■  v'»2  cross  section 
is  2.5a  .  See  the  caption  to  Figure  2  for  further  details. 
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I.  Introduction 

The  quenching  of  electronically  excited  Na  by  Hj  is  a  prototype  for 
electronic-to-vibrational  energy  transfer  in  molecular  collisions.  This 
process  has  been  widely  studied  both  experimentally^  ''  and 
theoretically .6  11  tjjeoret£cai  interpretation  of  the  quenching  process 

2  2 

has  focused  upon  the  states  of  A'  symmetry;  the  lowest  8^  state  has  a 

well  and  exhibits  an  avoided  crossing  with  the  ground  state  surface. 

Previous  dynamical  studies  indicate  that  these  two  electronic  states  can 

qualitatively  describe  the  quenching  process. A  quantitative 

description  of  this  process  requires  accurate  ab  initio  potential  energy 

surfaces  for  the  electronic  states  of  importance  and  the  nonadiabatic 

coupling  terms  between  these  electronic  states.  Ab  initio  calculations  of 
9  7—0 

the  Ak'  states  have  been  made;  however,  the  nonadiabatic  coupling  terms 
have  not  been  reported. 

Recently  a  diatomics-in-molecules  (DIM)  potential  energy  surface  has 
been  presented  for  the  three  lowest-energy  potential  energy  surfaces  of 
NaH^.^*^  The  two  lowest-energy  surfaces  agree  well  with  the  work  of 
Botschwina  and  Meyer.  ’  A  major  advantage  to  using  the  DIM  formalism  is 
that  it  provides  a  global  representation  of  the  potential  energy  surfaces 
and  the  couplings  between  the  surfaces. 

There  has  been  a  large  amount  of  work  on  developing  methods  for  the 
quantum  mechanical  treatment  of  electronic  transitions  in  atom-molecule 
collisions.^  ^  The  early  work  of  Zimmermann  and  George^  compared  the 
use  of  diabatic  and  adiabatic  representations.  In  the  adiabatic 
representation  coupling  between  different  electronic  states  introduces 
first-derivative  terms  in  the  coupled-channel  equations.  Zimmermann  and 
George  solved  the  adiabatic  coupled-channel  equations  by  transforming  the 
N  second-order  differential  equations  into  a  set  of  2N  first-order 
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differential  equations.  Baers  approach  to  the  problem  sought  to  obtain 
coupled  channel  equations  that  could  be  solved  using  the  efficient 
computational  algorithms  developed  for  second-order  differential 
equations.  In  his  approach  a  transformation  is  made  from  the  adiabatic 
representation  to  a  diabatic  representation  in  which  all  coupling  arises 
through  off  diagonal  elements  of  the  potential  matrix.  This  approach  has 
been  applied  to  collinear  atom-diatomic  collisions  but  is  limited  to 
systems  in  which  only  two  electronic  states  are  important.  More  recently 
Rebentrost  and  Lester^  have  used  a  similar  diabatic  transformation  to 
treat  the  nonreactive  F  +  H_  problem  in  three  dimensions. 

r 


The  previous  quantum  mechanical  calculation  on  Na  +  1^  by  McGuire  and 

Bellum6  used  the  two  lowest-energy  adiabatic  potential  energy  surfaces  as 

an  approximation  to  the  diagonal  diabatic  potential  energy  surfaces.  The 

off-diagonal  terms  of  the  diabatic  potential  were  modelled  as  gaussian 

terms  centered  on  the  location  of  the  minimum  of  the  energy  difference  on 

the  crossing  seam.  Cross  sections  were  caluclated  within  the  Infinite- 

1 8 

Order-Sudden  (IOS)  approximation;  however,  the  integral  over  the  IOS 
angle  was  approximated  using  only  the  90  degree  approach  of  the  Na  to  the 
molecule. 

In  the  present  calculation  we  describe  a  new  method  for  calculating 
electronic  transition  cross  sections  using  adiabatic  potential  energy 
surfaces  and  the  nonadiabatic  coupling  terms  as  input.  In  this  approach 
we  use  a  mixed  adiabat ic-diabatic  representation  in  which  the  motion  in 
the  diatomic  internuclear  distance  is  treated  diabatically  and  the  motion 
in  the  Na  to  center-of-mass  of  distance  is  treated  adiabatically .  We 
use  the  IOS  approximation  and  converge  the  calculations  with  respect  to 
the  numerical  integral  over  IOS  angles. 
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In  section  II  we  present  the  details  of  the  new  theoretical  method 
and  in  section  III  we  present  details  of  the  computational  procedure. 
Section  IV  contains  results  of  the  calculations  and  section  V  presents  a 
discussion. 


II .  Theory 

We  consider  the  collision  of  an  atom  A  with  a  diatomic  molecule  BC, 

where  R  is  the  distance  from  A  to  the  center  of  mass  of  BC,  r  is  the  BC 

internuclear  distance,  and  y  is  the  angle  between  the  R  and  r  vectors. 
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Within  the  IOS  approximation  the  total  Hamiltonian  is  given  by 


|.T  t 

'IOS 


h2  a2  p 
“  m]R  W2- 


h2£,(2+l) 


int 


(1) 


where  y  is  the  reduced  mass  for  the  R  motion.  The  Hamiltonian  for  the 
internal  motion  is  given  by 


int 


hi  32  • 


+ 


h2j ( j+l ) 

^BCr 


+ 


(2) 


where  y gj.  is  the  reduced  mass  for  the  r  motion.  The  Born-Oppenheimer 

electronically  adiabatic  wavefunctions  are  eigenfunctions  of  the 

electronic  Hamiltonian  H  . 

el 

Hel  <t>nA^-’r,R,Y^  =  V^(r,R.y)<ti^(£, r,R,y)  (3) 


where  x  is  the  collection  of  electronic  coordinates  and  Vna^(r,R,y)  is  Che 
electronically  adiabatic  potential  energy  surface  for  electronic  state  n. 
Within  the  IOS  formalism  the  internal  and  orbital  angular  momenta,  j  and 
are  constants  and  are  not  coupled  with  the  electronic  angular  momenta. 
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The  total  wavefunction  for  an  initial  state  ctg  is  expanded  in 
eigenfunctions  of  the  internal  Hamiltonian 


(x,r,R,y)  =  R_1£^(x,rtR,Y)  J  (R,Y> 

o  a  o 


where 


Hint  ^(x.r,R,y) 


(4) 


(5) 


The  wavefunctions  for  the  internal  degrees  of  freedom  (x  and  r)  are 
obtained  as  follows.  For  a  fixed  R  and  y  a  transformation  is  made  to  a 
representation  which  is  P-diabatic^  with  respect  to  the  r  motion 


id  A. 


aA, 


(6) 


This  transformation  is  defined  by 
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3r  unn 


(r,R,y) 


L  un(,nl(r,R,y) 


(7) 


where 


(S) 


and  the  subscript  on  the  matrix  element  denotes  the  variables  integrated 
over.  The  internal  wavefunctions  'l»a(x,r >R»Y)  are  expanded  in  this  P- 
diabatic  basis 


^(x_,r,R,Y)  = 


r_1  E  <tvfW,R,Y)  PV(D  C;Jv  a(R,Y) 
nv  ’ 


(9) 


where  P  (r)  is  a  vibrational  basis  function.  Substituting  (9)  into  (5) 

* 

and  closing  from  the  left  with  /  dx  f  dr  r2  r"1  <J>dA  UtR,y)  p*(r) 
yields  the  eigenvalue  equation 

H^(R,y)  C^(R,y)  =  C^R.y)  E^(R,y)  q 


where  the  matrix  elements  of  the  internal  Hamiltonian  are  given  bv 


Hnv.n'v'(R» 


Y>  =  /dr  p*(r )  j  f-  JlL  1  * 

J  v  1  L  2ubc  aF7  IK r2  J  6nn- 


+  C<r. 

dA  •  . 


R,y)  |  p^Cr) 


and  the  potential  matrix  V  ,  is  given  by 


VdA(r,R,y)  =  ijT(r,R,y)  VaA(r,R,y)  u(r,R,y) 


The  close  coupling  equations  for  a  fixed  y  are  given  by 

r  r  * 

J  dx  Jdrr2  ^((x,r,R,y)  [h}Js  -  E  ]  4^(x,r,R,y)  =  0 

Substituting  equations  (1)  and  (4)  into  (13)  gives 

E!  r.  hi  1L  ^  h2Z(Z+l)  .  1 

|  L  ^U  JR*  +  — +  Ea(R,y)  .  Ej  6a,a 

~  T  Fi’a(R^}  -  STCa'afR»Y>  j  xiy(R.Y)  = 


The  coupling  terms  are  defined  by 


Fa'afR’ y)  *  <  *J>(».r.R.T*  I  |jj  I  *J(x.r.R,y)  ^  ^ 

Ca'c.(R’v)  =  <  ^,(x,r,R,T>  |  1^  |  tJ(x,r,R,y)  > 


h 


Using  equation  (9)  these  can  be  rewritten  as 


Fa'c.<R-Y>  *  Eiciv,o,',R'V)[^„.<R.V)^,  (R.y) 

nvn' 

+  6nn'  Iff  Cn‘v,a(R,Y)] 


where 


Fnn,<R’Y>  *  <  C'^.y)  I  Iff  I  «V<i.R.Y>  > 


<Va<R’Y>  = 


£nrfv,C<R'Y)  [  Cnn'<R’Y)  lVv,a(R'Y> 

*  2  Fnn'  Iff  C^'v,a(R,Y>  +  5nn'  Iff7  Cn'v,(l<R’Y,]  <19) 


where 


C'(r-y)  ■  <  ^A<i’R’Y)  I  IffT  I  C  >, 


The  radial  wavefunctions  xta  (R,Y)  are  subject  to  the  standard 

0 
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boundary  conditions 


xaa  *R,y*  ^=0  "  0 


X““o<R’Y)  rU"  kj*  {  5a%  e*p[-HV  - 

■  €  <y)  «p[‘<v  -  *"«>]} 


where  the  asycptotic  wavenutr.ber  is  defined  by 


k,a  »  £  tE  -  Ei<R-Y)]  'i 
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Opacity  functions  can  be  defined  for  each  value  of  j  and  Z 


=  [  dy  siny  I  <$ 

a0tn  j 


-  S^2,  (y) 

aaQ  aa0 


(24) 


and  the  total  cross  section  for  transition  from  initial  state  j,  to  final 

IQ 

state  a  summed  over  all  final  i's  is  7 

o 


~~  ?  p““o 


(25) 
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III.  Calculational  details 
A.  Vibrational  basis  functions 

The  vibrational  basis  functions  are  chosen  to  be  eigenfunctions  of 
the  asymptotic  diatomic  Hamiltonian.  For  the  Na  +  system,  all  of  the 
electronic  states  considered  correlate  asymptotically  with  Hj  in  its 
ground  electronic  state.  Within  the  DIM  formalism,  all  electronic 
coupling  vanishes  for  Na  infinitely  separated  from  and  therefore  the 
asymptotic  diatomic  Hamiltonian  is  simply  that  for  a  vibrating,  rotating 
diatomic  molecule.  The  vibrational  basis  function  are  defined  by 


+  h2.T(j+l> 
2uBCr 


+ 


p  (r) 
Kv 


p  e  (r) 
v  v 


(26) 


where  Vg^Cr)  is  the  asymptotic  BC  vibrational  potential.  The  basis 
functions  are  expanded  in  a  harmonic  oscillator  basis,  (h^Cr)}, 


Pv(f) 


E 

k 


hk(r) 


kv 


(27) 


All  integrals  over  r  are  done  by  Gauss-Kerraite  quadrature;  both  for  the 
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integrals  necessary  to  define  Py(r)  and  the  integrals  necessary  in 
computing  'l>a(x,r,R,Y) .  The  numerical  values  of  the  vibrational  basis 
functions  are  stored  on  the  grid  of  quadrature  points  (r^j^to  be  used  in 
subsequent  numerical  integrals. 

B.  P-diabatic  transformation 

Equation  (7)  for  the  transformation  matrix  u  is  solved  by  the  Magnus 
method.^  The  transformation  matrix  is  only  needed  at  the  Gaussian 
quadrature  points  and  we  approximate  it  as  follows 


— (rm+l,R’^1  = 


exp 


-(r 


m+1 


-  r  ) 
rn 


u(rm,R,y) 


(25) 


The  choice  of  u  at  the  first  grid  point  is  arbitrary  and  the  final  results 
of  the  calculations  are  independent  of  this  choice.  The  exponential  of  a 
matrix  is  evaluated  by  the  method  of  ref.  20. 

C.  R-matrix  propagation 

Equation  (14)  is  not  solved  directly  for  the  radial  wavefunctions , 
instead  we  use  the  R-matrix  propagation  method  to  obtain  a  global  R 
matrix  for  relating  the  radial  wavefunction  and  its  derivative  a  large  R 
values.  The  S  matrix  is  then  obtained  from  the  R  matrix.  The  method  for 
treating  electronic  transitions  within  the  R-matrix  propagation  method 
when  the  input  is  in  an  adiabatic  representation  has  been  previously 
presented.^  In  brief  review,  the  internal  basis  functions  .ijja(x,r ,R,y)  are 
chosen  to  be  independent  of  the  radial  coordinate  R  within  each  sector. 
Propagation  across  a  sector  can  be  easily  expressed  in  terms  of  the 
eigenvalues  of  the  internal  Hamiltonian,  E'^(R,y).  The  coupling  between  the 
internal  states  arises  at  the  sector  boundaries  in  transforming  from  the 
internal  basis  functions  in  one  sector  to  those  in  the  adjacent  sector. 

It  is  the  defintion  of  this  sector  transformation  matrix  which  is  the  most 


-i  -  ■  -  >-■■■»—  ^  ^  «-!*-■  a—  >- 


r \r 
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difficult  aspect  of  applying  the  R-matrix  propagation  method. 

The  transformation  matrix  from  sector  i  to  sector  i+1  is  defined  by 


T^+1(Y)  =  <  ^(x,r,Rj,y)  Uj,(x,r,Rj+1,Y)  >i>r 


where  R^,1  is  the  center  of  sector  i.  It  has  been  shovm  that  this  can  be 


approximated  by 


Tia'+1(Y)  =  exp 


wj+1  -  r‘>  *  r*"*),y] 


where  F(R,y)  is  defined  in  eq.  (17).  We  use  an  alternate  method  to  obtain 
T*’*+*(Y).  Substituting  eq.  (9)  into  eq.(30) 


t^+1(y)  =  2  cl 


,a<RC,Y)  JdT  pt^r)Jd-  *nA  (*>r»Rc»Y) 


<J>n»(x,r,R^  1,y)  Py.(r)  Cn'v',a'^RC  ,Y^ 

dA 

However,  the  P-diabatic  basis  functions  <t>n  (x,R,y)  are,  by  definition, 
independent  of  r.  Using  this  fact  and  substituting  eq.  (6)  we  obtain 


4'*1(y)  *  X),  CJv,c(RJ'F>  u»n(r'RCY) 

nun  •  mm ' 


Since  the  l.h.s.  of  eq.  (32)  is  independent  of  r  we  can  choose  any  value 
for  r  in  the  r.h.s.  of  this  equation.  We  use  eq.  (32)  where  we 
approximate  the  overlap  of  the  adiabatic  electronic  state  between  sectors 


=  <  <fe*A(x.*r>Rc«y)  I  ,r,Rc+A’ Y^  >x 


aA,  _  Qi  +  1 


(33) 
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1.  Potential  energy  contours  for  the  three  lowest  adiabatic  potential 
energy  surfaces  V^^Cr.R.fc)  of  NaH2;  parts  (a),  (b),  and  (c)  are  for 
n  =  1,  2,  and  3,  respectively.  The  energy  contours  are  plotted  as  a 
function  of  the  H2  internuclear  distance  r  and  the  distance  from  Na 
to  the  center  of  mass  of  Hj  R,  for  a  fixed  angle,  =  90  degrees. 

The  contours  are  evenly  spaced  at  0.2eV. 

2.  Adiabatic  potential  energy  curves  VnaA(r,R,<)  and  diabatic  potential 

dA 

energy  curves  V  (r,R,&)  for  NaH,, .  The  curves  are  plotted  as  a 
function  of  the  H2  distance  r  for  fixed  Na-to-H2  distance  R  and  for 
fixed  angle  X  *  75  degrees.  The  diabatic  potential  curves  are 
defined  by  eq.  (12)  and  are  diabatic  only  in  the  r  motion.  For  parts 
(a)  and  (d)  R  =  3sq,  for  parts  (b)  and  (e)  R  =  4aQ ,  and  for  parts  (c) 
and  (f)  R  *  5ag .  Parts  (a)  -  (c)  show  the  three  lowest  adiabatic 
curves  (solid  lines)  and  the  three  diagonal  diabatic  potential  curves 
(dashed  lines).  Parts  (d)  -  (f)  show  the  off-diagonal  diabatic 
potential  coupling  curves:  1,2  (solid);  1,3  (short  dashed);  and  2,3 
(long  dashed). 
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1.  Potential  energy  contours  for  the  three  lowest  adiabatic  potential 
energy  surfaces  V^^Cr.R.fc)  of  Nal^;  parts  (a),  (b),  and  (c)  are  for 
n  *  1,  2,  and  3,  respectively.  The  energy  contours  are  plotted  as  a 
function  of  the  Hj  internuclear  distance  r  and  the  distance  from  Ha 
to  the  center  of  mass  of  Hj  R,  for  a  fixed  angle,  i  *  90  degrees. 

The  contours  are  evenly  spaced  at  0.2eV. 

2.  Adiabatic  potential  energy  curves  V  ^(r.R.O  and  diabatic  potential 
energy  curves  V^^Cr.R.'i)  for  Nal^.  The  curves  are  plotted  as  a 
function  of  the  Hj  distance  r  for  fixed  Ha-to-Hj  distance  R  and  for 
fixed  angle  X  ■  75  degrees.  The  diabatic  potential  curves  are 
defined  by  eq.  (12)  and  are  diabatic  only  in  the  r  motion.  For  parts 
(a)  and  (d)  R  *  3sq,  for  parts  (b)  and  (e)  R  *  4sq,  and  for  parts  (c) 
and  (f)  R  ■  5ag.  Parts  (a)  -  (c)  show  the  three  lowest  adiabatic 
curves  (solid  lines)  and  the  three  diagonal  diabatic  potential  curves 
(dashed  lines).  Parts  (d)  -  (f)  show  the  off-diagonal  diabatic 
potential  coupling  curves:  1,2  (solid);  1,3  (short  dashed);  and  2,3 
(long  dashed). 
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Figure  Kb) 


adiabatic  and  diagonal  diabatic  potential  curves 


Figure  2(a) 


diagonal  diabatic  potential  curves 


Figure  2(e) 


Figure  2(f) 


AN  ITERATIVE  APPROACH  TO  REACTIVE  SCATTERING 


This  appendix  presents  the  results  of  research  to  extend  the  iterative- 

•  1  2 

variational  method  ’  to  reactive  atom-molecule  collisions.  The  rationale  for 
this  work  is  the  following: 

•  The  standard  reactive  scattering  methods  limited  in  the  number  of  coupled 
equations  they  can  handle  (<200),  thus  limiting  accurate  close-coupling 
treatments  to  nearly  thermoneutral  reactions  involving  hydrogen 
abstraction. 

•  The  standard  reaction  coordinate  methods  are  limited  by  the  use  of  mass- 
weighted  coordinates,  and  hence  are  difficult  to  apply  to  small  skew- 
angle  systems.  This  comment  does  not  apply  to  so-called  hyperspherical 
coordinate  schemes,  although  the  limitation  to  systems  of  the  order  of 
100-200  channels  does. 

•  The  standard  close-coupling  theory  (using  Cartesian  coordinates)  contains 
nonlocal  potentials  to  account  for  rearrangement.  These  nonlocal 
potentials  arise  in  the  same  manner  regardless  of  whether  the  exchange 
involves  electrons  or  nuclei. 

•  Iterative  methods  have  been  very  successful  in  treating  nonlocal 
potentials  that  arise  in  electron  scattering  theory.  A  similar  iterative 
method  has  also  been  successful  in  solving  large  (>500)  sets  of  coupled 
equations  for  inelastic  atom-molecule  scattering. 

These  observations  suggest  that  an  iterative  method  might  form  the  basis 
for  a  theoretical  approach  to  reactive  scattering  that  can  be  applied  to  a  much 
wider  class  of  problems  than  can  refinements  of  current  reaction-coordinate 
methods.  The  rest  of  this  section  presents  a  discussion  of  our  preliminary 
research  in  this  area. 


B.l  INELASTIC  SCATTERING 

This  section  describes  unpublished  extensions  to  the  iterative  method  as 
it  applies  to  inelastic  scattering.  Consider  an  atom  A  colliding  with  a  target 
molecule  BC,  with  the  Hamiltonian  (refer  to  Figure  1) 


H  = 


.  Hi  72 

2u  VR 


Si 

2m 


V2 

r 


+ 


V(R,r,y) 


(1) 


where  Che  coordinates  r  and  R  correspond  to  the  BC  internuclear  vector  and  the 

vector  from  A  to  the  BC  center-of-mass ,  respectively,  y  is  the  angle  subtended 

by  R  and  r .  The  reduced  masses  are  given  by 

mB  mC  lo, 

m  =  -  (2) 

mB+  C 


=  mA  mB  mC 

H  is  more  convenient  when  the  potential  is  partitioned  into  a  term 
corresponding  to  the  BC  molecule  and  an  interaction  term 

Vg.  =  Lira  v(R,r,y)  (4) 

R-k» 


V  =  V  -  J 
int  BC 


so  that  the  moleculear  and  total  Hamiltonians  are 


n2 


25  Vr  +  VBC 


H  -  ££  1  3  p  +  M  „ 

“  -  -  2u  R  3R!  R  +  ^7  *  «B0  f  Vi„t 


where  L(R)  is  the  orbitel  angular  momentum  operator.  We  now  seek  solutions 
for  a  specific  total  angular  momentum  of  the  form 


!*f (J ,m)  = 


UJf(R) 

J.  I 


where  i_{n,j*S.)  and  f  ={n' ,  j  '  1 } .  y^  are  chosen  to  be  eigenfunctions  of  H, 


and  of  total  and  orbital  angular  Tiomentum. 


& 


tr- 


A  variational  principle  is  used  t  >  ~btiin  the  coupled  equations 
I  =  /drdR  ¥f*(H-E)  »f 

-  r*  e-  s  + /«  E°if  wij<B) 

J  i  J  ij 

’/«  £  I  "If  ♦/«*  £  »‘f  V*'  "if  -  £  "If  "If  c 

"J  1  •*  1J  ' 


(9a) 


(9b) 


(9«) 


wat;  re 


iisi 

S: 


r 

r-. 


r* 


_«2  (-jlU.+U 

W..(R)  =  f  - 

1J  |  2 


■S[ 

/• 


+  (e 


f  - «] 


+  I  dr  dR  ^int 


(10) 


The  best  solution,  the  one  that  leaves  I  stationary  with  respect  to  arbitrary 
variations  in  (endpoints  held  fixed),  is  given  by  the  solutions  to  the 
differential  equations 


uu-(R)  -^vr)  VR) 


(11) 


which  can  be  solved  by  a  large  variety  of  numerical  methods,  provided  the 
number  of  coupled  equations  is  limited  to  less  than  approximately  200.  We 
proceed  now  to  cast  these  equations  into  a  set  of  algebraic  equations,  to 
dramatically  increase  the  number  of  equations  that  can  be  simultaneously 
solved.  We  introduce  a  basis  for  the  radial  functions,  i.e. 


N 

VR>  *  £  *la<R)  Cio 

ct=1 


(12) 


The  variation  principle  then  leads  to  a  set  of  algebraic  equations  for  the 


4 


C 

l 


coefficients  C^a .  By  virtue  of  equation  (9c)  the  basis  functions  <|>^  need  only 
belong  to  the  class  of  continuous  functions  and  need  not  have  continuous  first 
derivatives;  i.e.  the  basis  may  be  peicevise  analytic.  This  last  observation 
suggests  that  the  finite  element  method  can  be  used  to  advantage  in  solving 
these  equations.  We  divide  the  range  of  R  into  a  grid  {  R1  <  r2  <...<  rn  }  and 
choose  the  <)>^a  such  that 


0  R  S  R 


a-1 


♦ia(*>  =  <  1  R  =  R 


(13) 


0  R  i  R 


a+1 


With  this  choice  =  ^f(Ra)«  Holding  the  endpoints  fixed,  the  variation 
principle  gives  the  Ncx(Ng-2)  algebraic  equations 


Nc  Ng 

ZZ 

j=l  8=1 


Kia,J8  Cj8  =  0 


2.3. 


N.,-1 


(14) 


where 

“ia.jB  ’J  dR  £  ho  *JB  SiJ  -  /®  Wij  »JB  <15) 

K  is  a  symmetric,  banded  supermatrix  (Figure  2)  because  of  our  definition 
of  the  d>^a  in  Eq.  (13).  Each  block  labeled  by  (a,S)  is  of  dimensions  NCXNC. 

Since  the  solution  must  be  regular  at  the  origin,  we  know  that  ■  o. 
For  any  choice  of  C.N  these  equations  may  then  be  solved. 

1  B 

It  is  important  to  stress  that  it  is  the  banded  symmetric  form  of  K  that 
makes  this  approach  attractive.  All  data  may  be  processed  sequentially,  thus 
avoiding  the  heavy  use  of  direct  access  I/O  that  characterises  the  previous 
implementation  of  the  iterative  method.  Therefore,  even  when  K  is  too  large  to 
fit  in  core,  the  solution  of  Eq.  (14)  is  of  comparable  effort  to  a  single 


numerical  integration  of  Eq.  (11). 

Additional  attractive  properties  of  this  approach  are  1)  there  are  no 
stability  problems,  even  for  solutions  beginning  with  incorrect  boundary 
conditions  at  ^^(Rh  )  and  2)  with  a  judicious  choice  of  ^a(R).  88  discussed 
below,  the  number  of  grid  points  can  be  considerably  reduced  over  the  number 
required  for  standard  numerical  integration.  Thus  this  approach  shares  a 
desirable  trait  of  the  more  popular  invar ient-embedding  methods  in  molecular 
scattering  theory*. 

The  C.„  are,  of  course,  not  known  in  advance.  However,  they  can  be 
iNg 

expressed  in  terms  of  the  others 

CiN_  s  Uif(RM  >  *  h  «if  “  elSif  (16) 

B  g 

where 


(21) 
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Substituting  Eq.(18)  and  Eq.(21)  into  Eq.(16)  leads  to 
Nc  JU 

BiNg,  j8  °  j  g  =  Xif(RNg) 

j=1  0=1 


where 


Bi»g,js  ■ ( t 1  *  *  ci  yi  f*. ft  wj ' j  ®j 


1SdR  +  SiNg.jB  <23> 


Eqs.(14)  and  (22)  together  result  in  a  set  of  NCX(K  "D  equations  for  Ncx(N  -1) 


unknowns 


K'  £  = 


where 

*38  =  5S«g  *jf<V  (25) 

and  Kf  has  the  structure  shown  in  Fig.  3. 

Eq . ( 24 )  can  in  principle  be  solved  directly  for  in  a  noniterative 

fashion.  However,  such  a  solution  is  unrealistic  for  large  N£  because  Eq.(22) 
destroys  the  banded  symnetric  nature  of  Therefore,  we  introduce  the 

following  iterative  procedure:  Let 

K'  *  K°  ♦  Kl  (26) 

where  JK°  is  banded  synnetric  and  ji’  is  aero  everywhere  except  in  the  last 


block  of  rows  (iN  ).  Let 
8 

K°  G1  -  % 


(27a) 


K°  2"  ♦  Kl  a"’1 


(27b) 


5°  -IdJ1 

i*1 

We  choose  the  d^  by  minimising 

I'  -  I K'  <5  n  -  el2 
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subject  tc  the  renormalization  constraint 

E  *1  • 1 

i 

This  leads  to 


where 


l  (j'1),. 

j  ~  iJ 


->i*  ♦j 

J  *  A  •  AJ 
J 


t1  =  Kl  (S1-  S1-1  ) 


Equations  (27)  through  (34)  are  iterated  until  convergence  is  achieved. 

An  excellent  choice  of  basis  functions  <J)^a  is  obtained  by  using 
eigenfunctions  of  the  equations 

4,ia(R)  =  Wii(R>  (35) 

These  equations  are  solved  and  the  integrals,  Eq.(15),  are  evaluated 
numerically  by  a  suitable  quadrature.  With  this  choice,  our  basis  functions 
always  have  the  correct  local  behavior,  thus  improving  convergence,  and  go  to 
the  proper  asymptotic  limit  when 

"ij  *  0  i*j  (36) 


This  swans  that  at  long  range  the  grid  spacings  may  be  ouch  larger  than  the 
local  wavelength,  and  that  fL.  may  be  kept  relatively  small,  even  in  the 
presence  of  coulomb  interactions. 
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B.2  REACTIVE  SCATTERING 

We  now  consider  an  extension  of  the  iterative  method  described  above  to 
three-atom  rearrangements  A  +  BC  -►  AB  +  C  but  exclude  breakup  A  +  BC  A  +B  +C 
since  we  are  still  limited  by  Lippmann-Schwinger  boundary  conditions.  We 


expand  our  solutions  as' 

N, 


ac 


(37) 


Y0f(J.M)  =  X^VRa)  — 

a=a,bc  i=1  Rot 

where  the  arrangements  atb,  and  c  are  shown  in  Fig.  4.  The  radial  functions 

have  the  boundary  conditions 


Ri®0OUai.8f(Ra)  =  eai  6if5afi  ‘  fai6aifi3f 
a 


Again,  we  expand  in  a  peicewiae  analytic  basis  *  4>aia» 
N 


(33) 


“al.Bf'V  *  £  V.11.1  °aio 
0=1 


Variation  of  the  C  ,  will  lead  to  the  set  of  equations 
at  a 


(39) 


X  K, 

Bjo 


aio.Bjp  C8jp  =  0 


o  =  2,3,  •••  N  -  1 

ag 


(40) 


| 
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We  use  integral  expressions  for  the  S-matrix  elements  to  provide  an  additional 

equation  for  each  arrangement  channel.  The  resulting  algebraic  equations  will 

be  solved  iteratively  as  in  Eqs.  (27)-(34). 

The  method  discussed  for  inelastic  scattering  in  the  previous  section  can 

be  taken  over  without  alteration  to  the  reactive  case.  There  are  differences 

in  application,  however,  that  will  now  be  considered. 

First,  the  matrix  elements  K  . _  0 will  now  form  a  supermatrix  of 

alo.BJp 

superaatrices,  the  block  structure  being  shown  in  Fig.  5.  The  diagonal  blocks 


with  a =8  contain  identical  integrals  to  the  inelastic  case  in  Eq.(lS).  The 
elements  of  the  off-diagonal  blocks  are  more  complex,  involving  overlaps  of  the 

A  ^  A 

nonorthogonal  functions  and  Xgj(rg»Rg)  *  These  blocks  are  not 
banded,  since  the  functions  $a^0(Ra)  will  in  general  have  nonzez'o  overlaps  with 
functions  for  other  arrangements.  The  off-diagonal  blocks  are  obviously  the 
exchange  interactions  and  are  matrix  elements  of  the  nonlocal  exchange 
kernals^.  The  associated  expressions  for  the  S-matrix  elements  are 
consequently  more  complex  than  Eqs.(18)-(21). 

In  addition,  the  matrix  K1  of  Eq.(26)  no  longer  consists  simply  of  the 
boundary  condition  Eq.(22).  Since  K1  is  no  longer  banded,  K°  can  be  chosen  by 
constructing  parallel  lines  to  the  diagonal  as  indicated  by  the  dashed  lines  in 
Fig.  S.  The  bandwidth  can  be  chosen  consistent  with  the  speed  and  size  of  the 
computer  being  used.  K°  will  include  everything  inside  the  band.  K1  will 
include  everything  outside  the  band  in  addition  to  the  rows  coming  from  the 
boundary  conditions  Eq.(22).  It  is  this  aspect  of  K1  which  simultaneously 
includes  the  boundary  value  problem  of  inelastic  scattering  and  the  nonlocal 
potential  problem  of  rearrangement  scattering  that  allows  its  treatment  within 
a  single  iterative  procedure.  Although  is  more  complicated  in  the  reactive 
case,  the  iterative  procedure  still  retains  its  sequential  processing  property. 
This  is  extremely  important  for  the  development  of  a  computationally  tractable 
method  on  the  new  generation  of  supercomputers,  with  their  enormous  fast-core 
memories  and  core-like  (fast)  online  storage,  and  vector  processors.  These 
extremely  large  and  fast  machines  will  enable  an  iterative  sequential-access 
procedure  such  as  discussed  here  to  handle  thousands  of  coupled  equations, 
which,  particularly  with  decoupling  techniques,  will  open  up  a  large  number  of 
processes  to  theoretical  investigation  that  arc  beyond  the  consideration  of  any 
existing  method. 
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FIGURE  CAPTIOUS 


1.  Coordinate  system  for  three-atom  nonreactive  scattering.  The  origin  of  R 
is  the  center-of-mass  of  the  BC  molecule. 

2.  Schematic  representation  of  the  banded  supermatrix  JK. 

3.  Schematic  representation  of  the  matrix  jK1  . 

4.  Coordinate  system  appropriate  for  rearrangement  scattering.  The  center- 
of-mass  of  each  diatomic  fragment  is  the  origin  of  the  corresponding 
vector  to  the  colliding  atom  in  that  channel. 

5.  Schematic  representation  of  the  supermatrix  occur ing  in  the  extension  of 
the  iterative  method  to  reactive  scattering.  Each  diagonal  block  is 
identical  to  the  corresponding  block  of  the  inelastic  case.  The  inelastic 
ft1  block  is  replaced  by  the  region  inside  the  heavy  dashed  diagonal  lines 
and  outside  the  banded  region,  and  is  chosen  as  appropriate  for  the 
problem  at  hand  and  the  available  computer  resources. 
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Figure  3 


